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Abstract 



Advances in computer technologies allow calculations in ever larger model spaces. To keep 
our understanding growing along with this growth in computational power, we consider a 
novel approach to the nuclear shell model. The one-dimensional harmonic oscillator in a box 
is used to introduce the concept of an oblique-basis shell-model theory. By implementing the 
Lanczos method for diagonalization of large matrices, and the Cholesky algorithm for solving 
generalized eigenvalue problems, the method is applied to nuclei. The mixed-symmetry basis 
combines traditional spherical shell-model states with SU(3) collective configurations. We 
test the vahdity of this mixed-symmetry scheme on ^^Mg and ^^Ti. Results for ^^Mg, obtained 
using the Wilthental USD intersection in a space that spans less than 10% of the full-space, 
reproduce the binding energy within 2% as well as an accurate reproduction of the low- 
energy spectrum and the structure of the states - 90% overlap with the exact eigenstates. 
In contrast, for an m-scheme calculation, one needs about 60% of the full space to obtain 
compatible results. Calculations for ^^Ti support the mixed-mode scheme although the pure 
SU(3) calculations with few irreps are not as good as the standard m-scheme calculations. 
The strong breaking of the SU(3) symmetry results in relatively small enhancements within 
the combined basis. However, an oblique-basis calculation in 50% of the full p/-shell space is 
as good as a usual m-scheme calculation in 80% of the space. Results for the lower p/-shell 
nuclei 44-48rp- '^^Cr, using the Kuo-Brown-3 interaction, show that SU(3) symmetry 
breaking in this region is driven by the single-particle spin-orbit splitting. In our study 
wo observe some interesting coherent structures, such as coherent mixing of basis states, 
quasi-perturbative behavior in the toy model, and enhanced B(E2) strengths close to the 
SU(3) limit even though SU(3) appears to be rather badly broken. The results suggest that 
a mixed-mode shell-model theory may be useful in situations where competing degrees of 
freedom dominate the dynamics, and full-space calculations are not feasible. 
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Chapter 1 
Introduction 



Selecting the right basis to perform calculations is a very essential step in analyzing any 
eigenvalue problem; it is especially true for many body quantum mechanical problems. When 
performing calculations, symmetries are also very important. Each of the fundamental quan- 
tities, such as energy {E ), linear momentum (p), and angular momentum (L), is conserved 
due to an exact symmetry of the physical space. It is well known that energy conserva- 
tion is due to time translational symmetry, linear momentum conservation is due to space 
translational symmetry, and angular momentum conservation is related to rotational sym- 
metry. Any mathematical description used in physics takes advantage of these symmetries 
and incorporates them explicitly. For example, in the Lagrangian (Hamiltonian) formalism, 
the Lagrangian (Hamiltonian) of the system is explicitly invariant with respect to the fun- 
damental symmetries, such as time translation, space translation and space rotation. Most 
examples of exactly solvable problems come from systems with some type of symmetry pQ. 

The notion of a stable equilibrium state, which is often related to an energy minimum, is 
another very important concept in physics. If Xq is an equilibrium point of the Hamiltonian 
function (if) for a classical particle, then H can conveniently be expressed in a Taylor series 
around Xq: 

H = + \k{x - xof + O (Ax") . 

This way, the harmonic oscillator described by the Hamilton H = + \kxP' turns out 
to be one of the most important model systems in physics with many applications [2]. 
SU {n) is the symmetry group of the n-dimensional harmonic oscillator, while Sp{2n, R) is 
the corresponding dynamical group. Thus, the SU{3) symmetry of the three-dimensional 
harmonic oscillator is a very important approximate symmetry of a system near equilibrium. 

Symmetries are very useful in the construction of shell-model structures in nuclear physics 
as well as in atomic physics. A shell-model structure is based on some exactly solvable limit 
of an effective interaction potential. An exactly solvable system allows for a well defined set 
of basis states. In particular, if bound states exist, they can be considered as single-particle 
levels of the system. Usually, a shell model assumes a mean field with which the particles 
of the system interact. For example, in atomic physics, the shell structure is mainly due to 
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the Coulomb field of the nucleus, while in nuclear physics the mean field is often taken to be 
the Hartree-Fock mean field. In this approach, the particle-particle interaction is assumed 
to be incorporated as much as possible in the average mean field. In particular, the nuclear 
spherical shell model is very successful in the description of nuclei [3] . Despite the enormous 
success of the spherical shell model, it is generally difficult to deal with nuclei in the middle 
of the shell (mid-shell nuclei) using this model. For such nuclei the collective degrees of 
freedom are very essential and the shell-model configuration space is very big. Therefore, for 
these nuclei, a shell model based on the collective degrees of freedom is more appropriate. 
Elliott's SU (3) model is useful for understanding the collectivity in light nuclei, up to A < 28 
(sd-shell) For heavier nuclei with A > 80, the pseudo-S'[/(3) version of Elliott's model 
is very successful in the description of the collective modes [S]. For these nuclei {A > 80,), 
the deformed Nilsson model is more accurate in the description of the single-particle levels 
than the simple spherical shell model [6]. 

At least in principle, collective phenomena, such as rotational spectra with strong B{E2) 
transitions, should be reproduced by the microscopic models. However, to do so using the 
spherical shell model, one needs sufficiently many particle configurations. Unfortunately, 
the dimensionality of the space grows combinatorially with the number of particles placed in 
the allocated levels. This binomial growth is a major computational problem. On the other 
hand, the SU (3) model allows for a good understanding of the collective nuclear properties 
in light and heavy mid-shell nuclei. However, for nuclei near closed shells, the spherical shell 
model is more favorable due to the dominance of the single-particle phenomena in these 
nuclei [7]. Therefore, it seems plausible to consider a hybrid-type calculation that uses these 
two models. In general, the two bases, the spherical shell- model basis and the SU{3) shell- 
model basis, will not be orthogonal to each other. Such a calculation can be considered as 
an "oblique" basis shell-model calculation [8]. 

The oblique-basis calculation for nuclei is the subject of the research presented here. 
Oblique-basis calculations are expected to be of a practical value in systems with competing 
degrees of freedom. For example, our study shows the relevance of the oblique calculation in 
the case of ^^Mg. For this nucleus, the single particle excitations described by the spherical 
shell model and the collective excitations described by the SU{3) shell model are important. 
When we combine the two bases, we obtain a significant gain in the convergence of the low- 
energy spectra towards the full space result. In particular, the addition of the leading-S'?7(3) 
irreducible representations (irreps) yields the right placement of the K = 2 band and the 
correct order for most of the low-lying levels. Indeed, an even more detailed analysis shows 
that the structure of the low-lying states is significantly improved through the addition of a 
few SU (3) irreps. 

The oblique-basis calculation will be an unnecessary numerical complication for systems 
where one of the excitation modes is dominant. For example, in the lower p/-shell nuclei ^^Ti 
and ■^^Cr, the spherical shell model gives a significant part of the low-energy wave functions 
within a few spherical shell-model configurations, while in the SU{3) shell-model basis one 
will need more than a few SU (3) irreps. This fact is mainly due to the strong breaking of the 
SU{3) in the lower p/-shell induced by the spin-orbit interaction In spite of the results 
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in the lower p/-shell, it is expected that in the mid-shell region some sort of SU{3) collective 
structure will gain importance If this is to happen, then the oblique-basis calculation will 
be an important alternative for calculating the structure of nuclei, such as ^^Fe and ^^Ni. 

Results of the shell-model calculations for lower p/-shell nuclei show that SU{3) symme- 
try breaking in this region is driven by the single-particle spin-orbit splitting. However, even 
though states of the yrast band exhibit SU{3) symmetry breaking, the results also show that 
the yrast band B{E2) values are insensitive to this fragmentation of the SU{3) symmetry; 
specifically, the quadrupole collectivity as measured by B{E2) transition strengths between 
low-lying members of the yrast band remain high even though SU (3) appears to be broken. 
Results for 44,46,48rpj ■^^Cr using the Kuo-Brown-3 two-body interaction [U] are given to 
illustrate these observations. 



^It was pointed by Chairul Bahri that the deformed Nilsson diagram for the p/-shell suggest a pseudo 
SU{3) symmetry. Another alternative could be a quasi-5'J7(3) symmetry. 
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Chapter 2 

The Nuclear Shell Model 



In some sense, the shell structure of nuclei is more complicated than the shell structure of 
atoms. The shell structure of atoms is due to the Coulomb force between the nucleus and the 
electrons. It may be a nice coincidence, but it is a fact that the Coulomb potential problem 
in quantum mechanics is an exactly solvable problem In the case of nuclei, the situation 
is more complicated. The reason is that there is no single source of a central potential. 
Instead, all nucleons are considered to act together, generating a mean field. Within this 
mean field, the problem is more tractable Here, we do not consider the problem of how 
to obtain the mean-field potential. Instead, we just use some general symmetry properties 
that a phenomenological potential and a realistic effective interaction should obey. These 
symmetry properties provide insight about the relevant single-particle basis within which 
one can consider the problem. 



2.1 Magic Numbers in Nuclei 

Maria G. Mayer's discussion of the magic numbers in nuclei has clearly demonstrated the 
nuclear shell structure associated with the independent-particle model for nuclei [TUj. In 
this model, each closed-shell configuration provides a convenient first approximation. In this 
approximation, one can assume that the system under consideration consists of a closed-shell 
core plus valence particles in a valence shell. This approach very successfully explains the 
ground state properties of nuclei [TT] . 

In order to understand and obtain qualitatively good results for the structure of the 
excited states, one has to consider a configuration mixing in the valence space. This usually 
leads to a very big model space. Therefore, a further truncation scheme is required. In this 
chapter, we will discuss two main approaches used in the nuclear shell model, namely the 
spherical shell-model truncation scheme and the SU (3) shell-model truncation schemes. 
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2.2 The Nuclear Interaction 



From a fundamental point of view, the problem of the relevant nucleon-nucleon interaction 
is very important. However, it is outside the scope of the research presented here. Even when 
one is provided with a good phenomenological nucleon-nucleon interaction, there is a lot of 
hard work to be done before one can finally set things up and calculate some experimentally 
meaningful results. Usually, a Hartree-Fock procedure is employed to reduce the many- 
particle Schrodinger equation to a single-particle Schrodinger equation with a self-consistent 
mean field. Once the single-particle states and energies are defined, then the n-particle 
configurations are formed using Slater determinants. Finally, a configuration mixing is used 
to take into account some of the residual interaction. This process may be simplified by 
using a phenomenological single-particle potential and a realistic interaction with a set of 
parameters adjusted to fit the experimental data. 

In this section, we consider a phenomenological interaction that contains some effective 
one-body and two-body potentials that are obtained from the original two-body nucleon- 
nucleon interaction: 

A A 
?=1 i^j s(i{valence particles} 

Ti is the kinetic energy of the i-th nucleon, V(|rj — rj\) is the two-body nucleon-nucleon 
interaction, ts is an effective one-body kinetic energy of the valence particles. Us is the 
effective mean-field potential, and Vres is the effective residual two-body interaction between 
the valence particles [3]. The effective one-body interaction H^^ = t + U provides a set of 
single-particle states: 

(x) = {t + U) 0, {x)=e4^ {x). 

The many-body wave function for a fermion system has to obey the Pauli principle. 
Thus, a fully antisymmetric combination, a Slater determinant, has to be constructed: 

01 (fi) 01 (fz) ■ ■ ■ 01 {Xn) 

02 (^l) 02(^2) ■■■ 02 (^n) 

0n (^l) 0n (^2) ■ ■ ■ 0n (^n) 

Here, the single-particle wave functions 0m (x^) correspond to the s-th particle in the m- 
th single-particle state with quantum numbers depending on the exact symmetries of the 
single-particle problem. Usually, these quantum numbers include angular momentum (j) 
and parity (vr). 

2.3 Hamiltonian in Second Quantized Form 

Given the single-particle levels, one can simplify the notation by going from the coordinate 
representation of the single-particle levels to an occupation representation. This process is 



* (fi, ....,f„) = det 
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often called a second quantization since the wave functions are constructed from appropriate 
creation/annihilation tensor operators acting on a vacuum state: 

Here, a stands for other quantum numbers, such as harmonic-oscillator shell numbers, spin 
and isospin labels. The vacuum state |0) is a reference state on which everything else is built. 
The vacuum state |0) may have a different meaning depending on the quantum labels of the 
annihilation operators. The annihilation operators usually define the vacuum as follows: 

flajm |0) =0. 

For example, if a^^-^ and a^^jra represent some real particles, such as fermions, then clearly 
the vacuum state is a state of no particles at all. If (i^j^ and a^jm represent the valence 
nucleons, then the vacuum state |0) would represent the closed-shell core. In the forthcoming 
chapters, we consider |0) to represent closed-shell nuclei. For example, ^^O is the closed-shell 
nucleus when we study nuclei in the valence sd-shell; ^°Ca is the closed-shell nucleus when 
we study nuclei in the valence p/-shell. 

In this second quantized form, the effective Hamiltonian is: 

H = ^Eiaftti + ^ Vij^kidtajakai. (2.1) 

Here, Ei are single-particle energies derived from excitation spectra of one valence particle 
system, i.e. ^^O in the case of the sd-shell. The Vki^ij are two-body matrix elements derived 
from an initial approximation, which are improved by a data fitting across the range of 
nuclei in consideration. For example, in the case of the sci-shell we would use the 63 two- 
body matrix elements obtained by Wildenthal [T2] . 



2.4 Spherical Shell Model for Nuclei 

We have already mentioned the independent-particle model [TT]. This model uses the 
harmonic-oscillator potential as an effective single-particle potential for nucleons [13] plus a 
spin-orbit interaction that provides for the correct shell closure [H]. In addition, there is a 
strong pairing part in the two-body interaction. The pairing interaction and the quadrupole- 
quadrupole interaction |Tlj are essential parts of the two-body interaction. 



2.4.1 Single-Particle Basis 

In computations based on the independent-particle basis, we use a phenomenological 
Hamiltonian (2.1) with single-particle levels labeled by the harmonic-oscillator quantum 
numbers nljm as follows: 

• n is the harmonic-oscillator shell. 
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• / is the angular momentum quantum number, 

• j = / ± I is the total angular momentum of the nucleon with spin 1/2, 

— * 

• m is the third projection of the total spin j. 

Within the above labeling scheme, the single-particle wave functions in the coordinate 
representation have the form: 

<Pnisjm (x) = {x\nlsjm) = ^ {Imi, sms\jm) R^i (r) Y"/^, (9, Lp) Xms- (2-2) 

m=mi+ms 

Here, {Imi, sms\jm) stand for the Clebsch-Gordan coefficients of SU (2), Rni{r) are the 
radial wave functions, Yjm; are the spherical harmonics, and Xms are the internal spin^ 

wave functions for nucleons. 



2.4.2 Many-Particle Basis 

In the occupation number representation. Slater determinant states are constructed from 
ni...nk nucleons by means of the fermion particle creation operators al : 



ni 



•••^fc)=nK)"|0)' (2.3) 



where the operators af and obey a Fermi algebra: 

aja^ + a^al = 0, 

a^aj + ajaj' = 6ij. 

Here, the labels of the operators and correspond to some specific quantum labels nlsjm 
of the spherical single-particle wave functions (2.2). 



2.4.3 Configuration Truncation and the M-scheme Basis 

Based on the independent-particle model, one can make an initial approximation to the 
wave functions of nuclei. This approximation uses the lowest energy configuration [ni_..., n^], 
where rij is the number of identical particles placed in the z-th orbital subject to the condition 
< n < 2j + 1. The energy of such a configuration is given by the expression ...,nu] = 
^^SiUi. It is immediately clear that in general there would be some degeneracy. Thus, the 
proper description of the excitation spectrum would need the two-body part of the interaction 
to lift this degeneracy. However, even then, using only the few lowest energy configurations 
is not sufficient to describe properly collective excitons in the mid-shell nuclei. 
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For heavy mid-shell nuclei, one needs to include a significant number of configurations. 
One way to proceed and include many configurations is to consider many-particle states 
with good J and Mj via SU{2) coupling within each configuration. Codes based on this 
approach usually rely heavily on 3j, 6j, and higher SU{2) symbols [I5l [16]. Since these j- 
symbols are calculated repeatedly, an efficient SU (2) package and a smart way to store often 
used coefficients are very essential. Recently, an SU{?>) code using the same strategy has 
been successfully developed [IT]. This code relies on a very efficient data storage technique 

An alternative computational method is the M-scheme approach [19]. In this approach, 
instead of using states with good J and Mj, one uses only states with good Mj and lets the 
Hamiltonian select the states of good J. Diagonalizing the Hamiltonian in such a basis results 
in a few of the lowest energy eigenstates. The M-scheme set of states is convenient since 
Mj is an additive quantum number. In order to provide for good total angular momentum 
(J), one has to include all states of fixed Mj within a given configuration. This method 
relies heavily on large matrix diagonalization algorithms. One such algorithm is the Lanczos 
algorithm which is very fast and efficient [20] • The Lanczos algorithm is a cornerstone of the 
modern M— scheme shell- model codes |21j . 

To illustrate the spherical shell-model truncation scheme, we consider ^^Mg. For this 
nucleus, the lowest configuration providing the initial approximation to the ground state is 
Qs^Qp^'^Qd^l^'^^s^i^Qd^l^. Here, Os^Op^^ is the core nucleus ^^O; the valence space is 0^5/2 lsi/2 
0(^3/2 with the lowest configuration of 8 particles, 4 protons + 4 neutrons, in the 0(^5/2 • If we 
explicitly write down a jj coupled state with good J and Mj within the Orf^yg configuration, 
then we would see that all the states with a fixed total Mj within the Od^^g configuration con- 
tribute to this state with good J and Mj. Since the Hamiltonian {H) respects the rotational 
symmetry, its eigenvectors must have good J and Mj values. Therefore, diagonalizing H in 
the space of all the states with fixed Mj within the Orf^^g configuration will automatically 
produce eigenstates with different J values and same Mj values. 

Usually, one has to include many configurations by using some selection principle. Often, 
the selection scheme uses the energy of the configurations. In this scheme, one includes only 
configurations that are within some range AE relative to the lowest energy configuration. 
Another selection scheme, which we use for the present study, considers the number of 
particles excited out of the lowest energy configuration into the full harmonic-oscillator 
shell. This selection scheme takes into account possible collective pair excitations when 
applied with two and four particle excitations outside of the lowest energy configuration]^ 

"'^Recently, it has been shown that one can successfully extrapolate some observables, such as energy 
eigenvalues, quadrupole moments, B(E2) transition strengths and Gamow- Teller transition strengths, using 
successively bigger truncation spaces. For more details see nucl-th/0203012 by Mizusaki and Imada and 
nucl-th/0112014 by Zelevinsky and Volya. 
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2.5 The SU(3) Shell Model for Nuclei 



If one considers a system near equilibrium, then it is possible to approximate its potential 
with a harmonic-oscillator potential. Since the symmetry group of the three-dimensional 
harmonic oscillator is SU{3), it is plausible to use SU{3) basis states. In this section we 
discuss the SU{3) shell model. We begin with a review of Elliott's SU{3) model [Ij. In 
particular, we present two single-particle labeling schemes, the spherical and cylindrical 
labeling scheme. Then, the structure of a general SU{3) irrep in the cylindrical labeling 
scheme is given. Next, we describe the SU{3) truncation scheme which is based on SU{3) 
invariant two-body interactions. We conclude the section with a brief discussion of the SU{3) 
breaking interactions. 

2.5.1 Labeling of the States in Elliott's SU(3) Model 

In this section we review group theoretical concepts that are important to the development 
of the theory and introduce SU{3) conventions adopted in our discussion. We consider the 
physical reduction, SU{3) D 5*0(3), and the canonical group reduction, SU{3) D f^(l) ® 
SU{2), with their respective labels. 

First we consider the physical group reduction SU{3) D 5*0(3). This reduction yields a 
convenient labeling scheme for the generators of SU{3) in terms of 50(3) tensor operators. 
The commutation relations for these SU{3) D 5*0(3) tensor operators are given in terms of 
ordinary 5*0(3) Clebsch-Gordan coefficients (CGC) {jm, j'm'\j"m") 

— V2{lm, lm'\lm + m')Lm+m', 

-VQ{2m,lm'\2m + m')Qm+m', (2.4) 
3"\/l0(2m, 2m'\lm + m)Lm+m'- 

Here, Lm are generators of the angular momentum and Qm is an algebraic quadrupole 
operator. 

Within this reduction scheme, states of an SU (3) irrep (A, fi) have the following labels: 

• (A,/i) - 5f/(3) irrep labels, 

• / - total orbital angular momentum, which corresponds to the second order Casimir 
operator of 50(3), 

• mi - projection of the angular momentum along the laboratory 2;-axis, 

• k - projection of the angular momentum in a body-fixed frame, which is related to 
multiple occurrences of 50(3) irreps with angular momentum / in the (A,/i) irrep. 

Unfortunately, this scheme has only one additive label, namely mi, and in addition, there 
are technical difficulties associated with handling the k label. 



mi Qm' 
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Figure 2.1: Three-dimensional view of the (A,yu) SU{3) irrep. 

The labehng scheme for our study is the canonical group reduction, SU{3) D U{1)^SU{2) 
|22j . In this scheme Qo is the U{1) generator and the SU{2) generators are proportional to 
Lq, Q+2, and Q-2 Under the action of the generators of these U{1) and SU{2) groups, 
the remaining four generators of SU{3) transform like two conjugate spin [|] SU{2) tensors 
with e = ±3 values for Qq. In this scheme, states of a given SU{3) irrep (A,yu) have the 
following labels: 

• (A,/i) - SU{3) irrep labels, 

• e - eigenvalue of the quadrupole moment (Qo), 

• nil - projection of the orbital angular momentum along the z-axis {Lq), 

• Up - related to the second order Casimir operator of SU (2), which for symmetric (A, 0) 
irreps is simply the number of oscillator quanta in the (x, y) plane. 
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This canonical reduction, SU{3) D U{1) ® SU{2), has two additive labels, e (Qo) and mi 
(Lq) and the allowed values of these labels for fixed SU{3) irrep (A,/i) are given by [23] : 



e = 2X + fi-3{p + q) (2.5) 
Up = fi+{p-q) 
mi = Up — 2m 



where 0<j9<A, 0<g</i, and < m < rip 



2.5.2 SU(3) Truncation Scheme 



It should be pointed out that the quadrupole operator Q used in (2.4) is actually an 
algebraic quadrupole operator 



with 6^ = — . However, the matrix elements of reduce to the matrix elements of the 
physical collective quadrupole operator Q'^ within a single harmonic-oscillator shell. 

i 

In general, the set of operators Q'^ and L are part of a Sp (6, R) Lie algebra. Within the 
Sp (6, R) model, the Q'^ operators connect same parity harmonic-oscillator neighbor shells 



The algebraic realization of the SU{3) model has the advantage that one can easily 
connect the important collective operators with the algebraically significant SU (3) operators. 
An example of a significant SU{3) operator is the second order Casimir operator of SU (3): 

Cf(^) = ^(3L2 + Q.Q). (2.6) 

By using the generators of SU{3) as labeled by the physical reduction SU{3) D 30(3), 
we can easily write a general algebraic SU{3) Hamiltonian: 

H = Hosc + xQ-Q + ■ L + aC^^^^^ + hL ■ Q ■ L + c{L ■ Q) ■ [Q ■ L) + dL ■ S (2.7) 

Here, Hose is the harmonic-oscillator Hamiltonian with single-particle energies e„ = /iw (n + |) . 
The strength of the quadrupole-quadrupole interaction is x- The 'bare' classical moment of 
inertia is when the effective moment of inertia will depend on a, 6, c and d. The parameter 
a is related to the third order Casimir operator cf^''^'' of SU{3). The strengths of the other 
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5*0 (3) invariant interactions, denoted by b and c, contain third and fourth order products 
of the SU{3) generators relevant to the multiphcity of the 5*0(3) irreps within the physical 
reduction 5[/(3) D SO {3) [25]. If 6 and care such that bL-Q ■ L + c{L-Q) ■ {Q ■ L) ^-fK'^ + L^ 
where 7 is the strength of the i^'-band splitting, then the collective states in the SU{3) D 
50(3) chain labeled by \N[f]{\fi)KLSJMj) would provide a basis in which H is diagonal. 



The main advantage of using an algebraic Hamiltonian, such as (2.7), is its SU (3) symme- 
try. Therefore, an SU (3) invariant Hamiltonian (H) does not connect states from different 
SU (3) irreps. Since the Q-Q interaction is proportional to the O2 of SU{3), it can be used as 
an essential SU{3) truncation scheme. This scheme prescribes 02-ordered importance of the 
SU (3) irreps. In this scheme, one selects SU{3) irreps (A, fi) with O2 = A^+/i^ + A/i+3(A-|-/i) 
values close to the biggest possible C2 value. The irrep with the biggest possible C2 value 
is called the leading SU{3) irrep. The leading irrep often corresponds to a total spin 5 = 
configuration. This way the leading irrep becomes also the dominant irrep for the low-lying 
energy states because the strength of the L ■ 5 interaction is usually expected to be small. 
This is due to the strong spin-pairing which tends to bring 5 = lower in energy. However, 
the one-body part of the ^ ■ /j ■ Sj interaction can cause significant deviation in the dominance 
of the leading irrep. 



Expressing the SU{3) Hamiltonian (2.7) in a second quantized form (2.1) gives: 



H = hij (n+^jJ2 '^tc^i + \ Y^ix {ij \Q-Q\kl) + ...)ata^akai. (2.8) 

i i,j 

Here, the labels i are shorthand notation for the single-particle labels in the SU{3) shell- 
model scheme, that is, i Tro{ri,0)Klsjmj in the SU{3) D 50(3) chain or respectively 
i — > TTo{ri,0)npemisms in the 5f/(3) D 5f/(2) chain. As usual, r = 1/2 is the isospin 
quantum number with tq = ±1/2 for protons/neutrons respectively, and (r/, 0) is the SU(3) 
irrep corresponding to a given harmonic-oscillator shell n [rj = n). The remaining labels 
were discussed in the previous section on the SU(3) shell model. 

2.5.3 Interactions that Break the SU(3) Symmetry 

Degenerate single-particle energies are an essential ingredient for good SU (3) symmetry; 



this is clear from our discussion on the general algebraic SU{3) Hamiltonian (2.7) and its 



second quantized form (2.8). However, we already discussed that the breaking of the single- 
particle degeneracy by the spin-orbit interaction is essential for the description of the correct 
nuclear shell closures in terms of the independent-particle model. Therefore, in case of 
a significant single-particle splitting, which is due to the orbit-orbit interaction ^ • If and 
the spin-orbit interaction li ■ Si, there would be a significant disturbance in the SU{3) 
truncation scheme. In this case, the spherical shell model described earlier would work and 
its truncation scheme could be used. 

Another SU(3) breaking factor is the pairing interaction. This interaction is a very 
essential short-range two-body nuclear interaction that can have significant impact on any 
5f/(3)-based calculations as well as on the spherical shell-model type calculations. Although 
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we have studied some effects of the pairing interaction in the sd-shell as well as in the pf- 
shell, we would rather not engage in this matter. We only mention that effects of the pairing 
in the context of the pseudo-5'[/(3) model have been studied before by C. Bahri [26], and 
currently we are considering incorporating the pairing effects within an oblique-basis type 
calculation via the broken pair model [3]. 
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Chapter 3 

Toy Model of a Two-Mode System 



The study of ^"^Mg, which will be discussed later in more detail, has successfully demon- 
strated the oblique-basis concept [8]. The quality of the results for ^^Mg are due to the near 
equal importance of the two basis sets used. On the one hand, the spherical shell-model 
basis is well-suited for description of the single-particle excitations; on the other hand, the 
SU (3) shell model puts an emphasis on the collective excitations in nuclei. These two modes 
are crucial for the ^^Mg example. 

In general, determining the relevant excitations is a cornerstone in the study of any 
system; in some sense this is the art of physics. Usually one basis works well for one system, 
but fails for another system. The reason is that in any general method, such as the variational 
method, perturbation theory, or fixed-basis matrix diagonalization, one needs to start with 
a good guess about the Hamiltonian and the states that describe the relevant excitation 
modes [27] . 

When applying perturbation theory, one is often concerned with a small perturbation of 
an exactly solvable limit of the full Hamiltonian [291 ED] ■ However, there are many examples 
when the relevant Hamiltonian has more than one exactly solvable limit [311 132]. This is 
a common situation when a dynamical symmetry group is used in the construction of the 
Hamiltonian [33], [31], [3S|- What shall we do if the system described by such a Hamiltonian 
is nowhere near any of the exact limits? In these situations, the problem may be better 
approached by using states associated with both limits. This set of states will form an 
oblique-mixed-mode-basis for the calculation. 

Taking into account the importance of the relevant energy scale of a problem and the wave 
function localization with respect to the range of the potential, the oblique-basis method 
can be taken beyond the idea of using two orthonormal basis sets. Specifically, one can 
consider a variationally-improved basis set starting with some initially guessed basis states. 
In the occupation number representation for the nuclear shell model, this variationally- 
improved basis method seems inapplicable!^ However, the method seems interesting because 

^In the occupation-number representation one assumes a fixed single-particle structure and then expands 
the states in the Slater determinants provided by this basis. From this point of view, there is no room for 
variationally-improved basis states since each Slater determinant is a single-integer machine word. 
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of its possible relevance to multi-shell ab-initio nuclear and atomic physics calculations. The 
method may also be related to some renormalization-type techniques. Therefore, a brief 
discussion of the variationally-improved basis and its possible applications is given in the 
Appendix. 

In this chapter some relevant mathematical notation and concepts used in the oblique- 
basis method are introduced. Specifically, we demonstrate the concept of the oblique basis 
on a simple two-mode system, the one-dimensional harmonic oscillator in a box. First, we 
discuss the concept and then the two exactly solvable limits of our toy model are briefly 
summarized. A qualitative discussion of the expected spectrum of the one- dimensional har- 
monic oscillator in a box is given. This is followed by an example spectrum and quantitative 
estimates. Some specific problems related to the structure of the Hilbert space will be ad- 
dressed. Finally, the main results will be discussed, especially a quasi-perturbative behavior 
and a coherent structure within the strong mixing region. 



3.1 Harmonic Oscillator in a One-Dimensional Box 

Let us start with an abstract two-mode system. For simplicity, we assume that the Hamil- 
tonian for the system under investigation has two exactly solvable limits, for example: 

H={l-X)Ho + XHi + \{l-X)H2. (3.1) 

Here, Hq and Hi are two exactly solvable Hamiltonians. This way, we have H ^ Hq m the 
limit A — 7- and H ^ Hi when A — ?■ 1. In the vicinity of these two limits we can approach 
the problem using standard perturbation theory. However, for A ~ ^ we have a very mixed 
system with unclear behavior which could be complicated further by an interaction H2 
between the natural modes of Hq and Hi. 



In the expression (3.1), A is introduced to simplify the discussion. In general, we have 
more than one parameter in the Hamiltonian. Often the exactly solvable limits are described 
as hypersurfaces in the full parameter space. It could even be that there are three or more 
exactly solvable limits. For example, the Interacting Boson Model (IBM) has three exactly 
solvable limits [2]. Another example with three exactly solvable limits is the commonly used 
nuclear schematic interaction. It has nondegenerate single-particle energies (£»), pairing 
{P~^P) two-body interaction, and quadrupole-quadrupole {Q ■ Q) two-body interaction: 

H = €iNi + GP+P -xQ-Q- 

Here, we consider the simplest two-mode system that is sufficiently close to the problem 
we have to solve for nuclei. The system under consideration consists of a one-dimensional 
harmonic oscillator in a one-dimensional box of size 2L l36i: 



H^l-r?^n(^^^e. (3.2) 
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where V^lq) is the confining potential which is zero for |g| < L and oo for |g| > L. This 
system has two exactly solvable limits. A more realistic model might consist of a three- 
dimensional harmonic oscillator and a square-well potential since these two potentials are 
known to be good starting points in the nuclear shell model [37] . 

The one- dimensional harmonic oscillator in a one-dimensional box model has been used 
as an example by Barton, Bray, and Mackane in their discussion on the effects of the distant 
boundaries on the energy levels of a one-dimensional quantum system |3H]. Also, some 
studies have already been done for the cylindrical symmetric system of a three-dimensional 
harmonic oscillator between two impenetrable walls [39j. However, the bi- modal structure 
of the problems has not been discussed in these studies. The essential two-mode regime 
of such problems has been studied in the context of a two-dimensional confinement of a 
particle in an external magnetic field by Rosas et al. |10]. Some authors have generalized 
the one-dimensional harmonic oscillator by introducing time dependent parameters in the 
Hamiltonian [4lJ and have recognized the two limiting cases of a free particle and harmonic 
oscillator. The infinite square well and the harmonic oscillator have been considered as the 
two limiting cases of a power-law potential within the context of wave packet collapses and 
revivals |l3l Il2]. Here, we focus our study on the bi-modal structure of the one-dimensional 
harmonic oscillator in a one-dimensional box. 



The first limit of the toy model (3.2) is = 0. This is a free particle in a one-dimensional 
box with size 2L : 

Ho = ^p' + VLiq). (3.3) 
The eigenvectors and energies are labeled by n = 0, 1, ... and are given by the expressions: 

y cos ((n + 1)^7-) if n is even 
*n(g) = <; , (3.4) 



^ sin (^(n + 1)||;) if n is odd 

^■^ = 2^ (<" + !) 2) (i 

This limit corresponds to extreme nuclear matter when the short range nuclear force produces 
an effective interaction well represented by a square- well potential [37]. We can think of 
this limit as a one- dimensional equivalent of a three-dimensional model where nucleons are 
confined within a finite volume of space representing the nucleus. 



The other exactly solvable limit of the toy model (3.2) is the harmonic oscillator in one 
dimension: 

1 o mcj^ 



In dimensionless coordinates 



1 ^ Iv ' p ^ pVmhw, 

mco 
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we have: 

H, = nuj\ {p' + ■ 

Thus the eigenvectors and energies are labeled by n = 0, 1, ... and are given by the expres- 
sions: 

En = hw (^n + ^ . 

Where Hn are the Hermite polynomials. This limit corresponds to the three-dimensional 
harmonic oscillator model for nuclei. 

In a one-dimensional toy model, the anharmonic oscillator with a quartic anharmonicity 
would be the appropriate counterpart of the Sp{6,R) shell model since the quadrupole- 
quadrupole interaction Q ■ Q goes as ~ and Q connects same parity harmonic oscillator 
shells. If we restrict the model space to only one harmonic oscillator shell, then we can use 
the algebraic quadrupole moment Q of Elliott [4] because within a single shell Q is the same 
as Q [2\. Thus for our study it is appropriate to consider the one-dimensional harmonic 
oscillator to correspond to the SU (3) shell model for nuclei. 



3.2 Spectral Structure at Different Energy Scales 

Often in physics the spectrum of a system is different for different energy scales. This 
usually reflects the existence of different excitation modes of the system. For the toy model 



Hamiltonian (3.2) we can clearly define three spectral types: 



Spectrum of a particle in a one-dimensional box (3.4) with quadratic dependence on n 
{En ~ n^). 



Spectrum of the one-dimensional harmonic oscillator (3.6) with linear dependence on 
n {En ~ n), 

Intermediate spectrum that is neither of the above two types. 



From Fig. |3.1| we expect that the particle in a box spectrum should be operative at high 
energies. These energies are energies where the box boundaries dominate over the harmonic 
oscillator potential. In this regime one can use standard perturbation theory to calculate 
the energy for a particle in a box perturbed by a harmonic oscillator potential. It can be 
shown that perturbation theory will give better results for higher energy levels. For n — )■ oo 
the first correction {SEj^) approaches the constant value of An estimate on when 

the perturbation calculations are feasible using — E^ » {n \ V\ n) gives: 

n » 2mVLV(3nV2). (3.7) 
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Figure 3.1: Two-mode toy system. The structure of the interaction potential of a particle 
in a one-dimensional box subject to a harmonic oscillator restoring force towards the center 
of the box. 



18 



30 



26 



22 



14 



10 



7f;7fl 38.5172 



^1,22^ _30.9981 



24.4686 



•1.8.-&1&7- 
14.3297 



3.0797 



ID HO 



Exact Perturbed 



E^=19J392 



AS. 



ID Box 



Figure 3.2: Spectral structure of the two-mode system ioi m = h = 2L/7i = 1 and u = 4. 



This analysis is confirmed by the numerical calculations shown in Fig. 3^ where the per- 
turbed particle in a box spectrum is really operative at n > 3 for the case of m = h = 
2L/7r = 1 and u = 4. 

The intermediate spectrum should be observed when the harmonic oscillator turning 
points coincide with the walls of the box. Therefore, the critical energy scale that separates 
the two extreme spectral structures is given by: 

= —L' (3.8) 

Notice that the constant energy shift in the energy of the high energy levels 5E\^y^ 

is one-third of the critical energy {Ec/3). 

At low energies, where the one-dimensional harmonic oscillator determines the classical 
turning points to be far from the boundaries, we expect to see the harmonic oscillator 



spectrum as shown in Fig. 3.2 The number of harmonic oscillator states that will be 



19 



observed is easily estimated using: 



J^c>K. ^^max-2~^^ 2 ^ 

It should be pointed out that there is a compatible number of levels, usually bigger than 
'^max' below the Ec corresponding to a free particle in a box: 

in in 2 muL? 
> E'f nZ. = ^ - 1. (3.10) 

However, these states are mixed by the harmonic oscillator potential toward the correspond- 
ing harmonic oscillator wave functions. 

Using the ratio of the ground state energies, 

Ef2/El^s. = 4mu;LV(n7r2), together with 



(3.9) and (3.10), the following spectral situations apply: 



For ^'^jI''^ > (f)^ there are levels below Ec corresponding to the harmonic oscillator 
and the free particle in a box such that E^^^ > E^^. However, only the harmonic 
oscillator levels are seen in the low energy spectrum. 



For (l) > '"'^^^ > I there are only the ground states E^^ and E^^ below Ec and 
El"" >E"o ' 



9-S- ' 9-s. 



• For I > '^"^^'^ > 1 there is only the ground state of the harmonic oscillator E^^ below 
Ec. 

Therefore, the smallest number of state^to illustrate the two mode spectra is the case of 



m = h = 1, L = 7r/2 and u = 4. With these parameters, formula (3.9) gives = 4. 5348. 
Thus one should see no more than 4 equidistant states as shown in Fig. 3.2 In Fig. 3.2 there 
are three clear equidistant energy levels that correspond to a harmonic oscillator spectrum. 

With respect to the critical energy Ec, there is a more explicit classification of the spectral 
structure: 

• Perturbed particle in a one-dimensional box spectrum for energies E » Ec such that 



(3.7) holds. 



One-dimensional harmonic oscillator spectrum (3.6) for energies Ec >> E such that 



(3.9) holds. 



Intermediate spectrum for energies E ^ Ec- 



^For simplicity we usually fix the parameters as follows: m — h = 1, L = 7r/2. 
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3.3 Toy Model Calculations and Results 



Despite the simplicity of the toy model (3.2), the harmonic oscillator in a box exhibits 



some of the essential characteristics of a more complex system. Our main interest is in 
problems associated with the use of fixed-basis calculations. In particular, one such problem 
is the slow convergence of the calculations [36] • If one can implement an exact arithmetic, 
one may not worry too much about the slow convergence when enough time, storage, and 
other resources are provided. However, numerical calculations are plagued with numerical 
errors that may grow significantly and render the results meaningless. From this point of 
view, a calculation that converges slowly may be compromised by accumulated numerical 
error. 



3.3.1 On the Hilbert Space of the Basis Wave Functions 

Before discussing the toy model using an oblique basis, it is instructive to discuss briefly 



the harmonic oscillator problem (3.5) using the wave functions for a free particle in a one- 



dimensional box ( 3.4[ ); and vice versa, solving the problem of a free particle in a one- 



dimensional box (3.3) using the wave functions for a particle in the harmonic oscillator 



potential (3.6) 



Due to the structure of the wave functions, there are some specific problems that need 
to be addressed. For example, using wave functions for a free particle in a one-dimensional 
box to solve the harmonic oscillator problem may not be appropriate especially for high 
energy states E » Ec- The problem is that any linear combination of wave functions 
with the same localized support, in our case the wave functions are localized within the 



box, will still be a function with the same localized support (see Fig. 3.3). That is, any 
linear combination of wave functions that are zero outside of the box is a function that is 
zero outside of the box too. Because the harmonic oscillator potential gets wider for higher 
and higher energies, any higher energy wave function must spread more than the previous 
one. Similarly, the spreading of the harmonic oscillator wave functions is responsible for the 
troubles that arise in solving the problem of a free particle in a one-dimensional box using 
the harmonic oscillator wave functions. The essence of these problems is in the structure of 
the corresponding Hilbert spaces. 

The influence of the boundary conditions on the properties of a quantum mechanical 
system has been recognized from the dawn of quantum mechanics. It is well known that 
some separable problems may re-couple due to the boundary conditions [15]. Some recent 
studies on the problem of confined one-dimensional systems using equations for relevant 
cut-off functions have been pioneered by Barton, Bray, and Mackane [38]. Their method 
has been further developed in a more general setting by Berman [H]. Other authors aim 
at variational procedures using simple cut-off functions [ITJ HH] or derive asymptotic esti- 
mates for multi-particle systems using the Kirkwood-Buckingham variational method [H]. 
Somewhat different approaches focus on shape-invariant potentials and use supersymmetric 
partner potentials to derive energy shifts and wave function approximations [35], as well 
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Figure 3.3: Spreading of the wave functions for the harmonic oscillator (blue) and particle 
in a box (red). 
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as sample-size dependence of the ground-state energy [50] • In the next few paragraphs we 
discuss the structure of the relevant Hilbert spaces when confinement is present. 



Harmonic Oscillator in the One-Dimensional Box Basis 



Now we consider the harmonic oscillator problem (3.5) using the wave functions for a free 



particle in a one- dimensional box (3.4). There are no difficulties for energies E « (3.8) 
where the harmonic oscillator potential is still within the box. However, for energies E » 
Ec the basis wave functions are localized only on the interval [—L,L]. Thus they cannot 



provide the necessary spread over the potential width (Fig. 3.3). This situation would be 



appropriate for the toy model (3.2) but not for the pure harmonic oscillator problem (3.5). 

One simple solution of the spreading problem is to continue the basis wave functions 
by periodicity. This way the necessary spread of the basis wave functions can be achieved 
and the new basis will stay orthogonal but must be re-normalized|^ However, these basis 
wave functions do not decay to zero in the classically forbidden zone. This means that 
some significant number of basis wave functions will be needed to account for the necessary 
behavior within the classically forbidden zone. 

Another alternative is to change the support domain corresponding to non-zero values 
of the function by stretching or squeezing it through a scaling of the argument of the basis 
wave functions, x — )■ xan/L. This way the support becomes [—L,L] [—«„,«,„]. Here, a„ 



is a scale factor for the n-th. basis wave function (3.4) estimated either from the width of 
the harmonic oscillator potentiaQ or determined by variational minimization. Either way, 
the new set of basis functions will be non-orthogonal. In general, there may be even a linear 
dependence. However, for the basis functions discussed here, linear dependence may not 
appear due to the different number of nodes for each wave function. The number of nodes 
(zeros) is not changed under the re-scaling procedure. While the potential width scaling is 
simpler, its applicability is more limited than the variationally-determined one. In general, 
the variational approach can be extended for much more general situations as discussed in 
the Appendix. 

• Particle in a Box in the Harmonic-Oscillator Basis 

Next, suppose we want to solve the problem of a free particle in a one-dimensional 



box [— L, L] (3.3) using the harmonic oscillator wave functions (3.6). The first thing to 
do is to change the inner product of the wave functions: (/, g) = /* (x) g (x) dx — > 

/-L /* (^) 9 (^) Then, it is immediately clear that the set of orthonormal harmonic os- 



cillator wave functions "^niQ) ( |3.6 ) will lose its ortho normality and even its linear indepen- 



^If one continues the wave functions to infinity, then there is a normalization problem. However, if the 
continuation is on a finite interval, then the functions can still be normalized. 

''The initial idea is to use basis states that have a spread compatible with the width of the potential in 
the energy region of interest, thus resolving the spectra only within that energy scale without calculating 
the lower energy states. Unfortunately, it does not seem to work since interference causes reduction of the 
wave spread and therefore drives the solutions towards the lowest eigenstate. 
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dencej^ However, this is not the actual trouble in such an approach]^ Neither the variational 
nor the potential-width wave function scaling will help to cure the loss of hermiticity of the 
physically significant differential operators, such as the momentum operator (p = —ih^) 
and the Hamiltonian operator {H = ^p"^)- This non- hermiticity is due to the behavior of 
the basis states at the boundary, mainly the non-vanishing of the wave functions at —L and 
L. For detailed analysis on the loss of hermiticity, we refer the reader to the Appendix. In 
order to recover the hermiticity of the differential operator it is sufficienlj^ to make sure 
that our basis wave functions vanish at the boundary points —L and L. For this purpose 
one can look at the nodes of each basis wave function and scale it so that its outer nodes are 
at the boundary points]^ Since the physical requirement that the wave functions have to be 
zero at the boundary is the cornerstone in quantizing the free particle in a one-dimensional 



box (3.4), it is not surprising that the nodally adjusted harmonic oscillator wave functions 



are very close to the exact wave functions for the free particle in a one-dimensional box as 



shown in Fig. 3.4 



In general, calculating the nodes of a function may become very complicated. To avoid 
problems with finding the roots, one can use the following techniqu^ the idea is to evaluate 
the value of the wave function at the boundary points, then shift the wave function by a 
constant to get zeros at the boundary, "^{q) — )• \E'(g) — \E' (L). This idea works well for 
even wave functions, but has to be generalized for odd wave functions by adding a linear 
term, — )■ "^{q) — (L) /L)q. Thus for a general function we can have: — )■ 



- (1 + q/L)"^ (L) /2 - (1 - q/L)-^ /2. In Fig. |3^ we have shown some of the 
resulting wave functions. Notice that this procedure gives a new wave function \1/ that is 
well behaved inside the interval [—L, L] and grows linearly with q outside the interval [—L, L]. 
This is in contrast to the behavior of the cut-off function /(g) obtained by Barton et al [38] . 
The function /(g) has L/q singularity at the origin (g = 0). The use of a cut-off function 
to enforce boundary conditions has been developed by Barton et al [3H] and Berman |13] 
and provides an interesting integral equation for the cut-off function. On the other hand, 
a simple cut-off function supplemented by a variational method seems to be very effective 

gziiiniiiH]. 

It should be pointed out that by using the above process one can set up and successfully 
run a modification of the usual Lanczos algorithm, to be discussed later, to solve for the few 
lowest eigenvectors of the free particle in a one-dimensional box through an arbitrarily chosen 
initial wave function. The major modification is to project every new function, = H^n-, 



■^The set of functions {q) with support domain restricted to [—L^L] and denoted by ^„(q; [— L,i]) 
may become linearly dependent if L is so small that there are more than one ^n{q\ [~L^L\) with the same 
number of nodes within [—L, L\. 

^The oblique basis type calculations described later can successfully remove the linearly dependent basis 
states in the process of handling the non-orthogonality of the basis. 

''Wave functions with the same value at ± L is the necessary condition; the wave functions should be zero 
only for an infinite potential at ± L. 

^From the nodal structure of the harmonic oscillator wave functions, given by the Hermite polynomials, 
it is clear that the first two wave functions (5*0 and 'Si) cannot be used since they have less than two nodes. 

^ This technique has been suggested by Professor A. R. P. Ran (private communications). 
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(a) (b) (c) 



Figure 3.4: Harmonic-oscillator trial wave functions adjusted with respect to the one- 
dimensional box problem: (a) adjusted according to the potential width E^^^"^ = =^ 
i^n = (1 + 2n), (b) nodally adjusted, (c) boundary adjusted using — > ^(q)—^ (L) (1+ 
q/L)/2-^{-L){l-q/L)/2 
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into the appropriate Hilbert space and subtract the components along any previous basis 
vectors. Only then should one attempt to evaluate the matrix elements of H related to 
the new basis vector that is clearly within the correct Hilbert space. This way, one has to 
double the number of scalar product operations compared to the usual algorithm where the 
matrix elements of H are calculated along with the complete re-orthogonalization of the 
basis vectors. 

3.3.2 Discussion of the Toy Model Results 



Having considered the main problems one may face in studying the simple toy model (3.2 ) 



we close the discussion with a sample spectrum for the case oi m = h = 2L/'k = 1 and 
w = 4. 



As one can see in Fig. 3.2 , the first three energy levels are really equally distant from one 



another and coincide with the harmonic oscillator levels as expected from (3.9). For these 
states, the wave functions are also the harmonic oscillator wave functions. The intermediate 
spectrum is almost missing. After the E^^ the spectrum is that of a free particle in a 
ID box perturbed by the harmonic oscillator potential. The oblique-basis type calculation 
reproduces the first eight low energy states within a 14-dimensional calculation, seven nodally 
adjusted harmonic oscillator states and seven states of a free particle in a box, while the 
fixed-basis calculation, using only the wave functions of a free particle in a one-dimensional 
box, requires 18 basis states. 

Due to the simplicity of the toy model, one does not find any big numerical advantage 
of the oblique-basis calculation compared to the calculations using the fixed basis of the ID 
box wave functions. There are two main reasons for this: (1) there is a sharp energy scale 
Ec that separates the two modes, (2) the spectrum above the energy Ec has a nice regular 
structure. 

The nice regular structure above the energy E^. results in a very favorable situation for 
the usual fixed-basis calculations since the dimension of the space needed to obtain the n- 
th eigenvalue grows as n + a. The parameter a is relatively small and does not change 
much in a particular region of interest. For example, the w = 16 calculations need only 
a = 15 extra basis vectors when calculating any of the eigenvectors up to the hundredth 
vector. The relatively constant value of a can be understood by considering the harmonic 
oscillator potential as an interaction that creates excitations out of the n-th unperturbed ID 
box state. Therefore, a is the number of ID box states with energies in the interval E^ and 
E^ + uP /2 |$„) where E^ is the n-th unperturbed ID box state energy. There is a 

fast de-coupling of the higher energy states from any finite excitation process that starts out 
of the n-th state. The fast de-coupling is due to the increasing energy spacing of the ID box 
spectrum. This results in a finite number of states mixed by the presence of the harmonic 
oscillator potential. Using the upper limit -Ec/3 on 5E\, one can easy estimate a: 
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Figure 3.5: Absolute deviations from the exact energy eigenvalues for u = 16, L = tt/2, h = 
m = 1 as a function of n. Blue circles represent deviation of the exact energy eigenvalue from 
the corresponding harmonic oscillator eigenvalue {AE = E^"''^^ — E^^), the red diamonds 
are the corresponding deviation from the energy spectrum of a particle in a ID box {AE = 
£jexact _ E^f^), and the green squares are the first-order perturbation theory results. 



^1 II) 

The sharp separation of the two modes allows for a safe use of the harmonic oscillator 
states without any rescaling. This is especially true when u is very big since then the low 
energy states are naturally localized within the box. Therefore, there is a clear shortcut: 
instead of diagonalizing the Hamiltonian in some ID box wave-function basis, one can just 
use the harmonic oscillator wave functions. 



Fig. 3.5 shows the absolute deviation {AE = i?^^"^* — i^estimate-j ^j^g exact energy 
spectrum for the case ofa; = 16, L = 7r/2, h = m = 1. Here, i^^'^**'"'**^ refers to the three 
energy estimates one cam make: the harmonic oscillator E^^, particle in a ID box E^^^, and 
the first order perturbation theory estimate considering the harmonic oscillator potential as 
a perturbation {E^^ + 00"^ /2 \^n))- There are about 19 states that match a harmonic 



oscillator spectrum which is consistent with the expected value from (3.9). After the n = 20 
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Figure 3.6: Relative deviations from the exact energy eigenvalues for u = 16, L = n/2, 
h = m = 1 as a. function of n. The blue circles represent deviation of the exact energy 
eigenvalue from the corresponding harmonic oscillator eigenvalue [AE/E = 1 — E^'^/E^^"-'^^), 
the red diamonds are the corresponding relative deviation from the energy spectrum of a 
particle in a ID box {AE/E = 1 — E^^ / E^"''^^) , and the green squares are the first-order 
perturbation theory results. 



level, the perturbation theory gives increasingly better results for the energy eigenvalues. 
Fig. 3.6 shows the relative deviation (1 — E<^stimate ^ j^exact^ ^^le exact energy spectrum for 
the case of a; = 16, L = 7r/2, h = m = 1. 

From these graphs, it seems that the transition region is somewhat absent since the 
first-order perturbation theory takes on immediately after the breakdown of the harmonic 
oscillator spectrum. Even though the first-order perturbation theory gives good estimates for 
the energy levels in this transition region, this is not a manifestation of a proper perturbation 
theory. Rather, it is a manifestation of a coherent behavior [51j- What actually happens in 
this region is a coherent mixing of ID box states by the harmonic oscillator potential in the 
sense of a quasi-symmetry discussed in the Appendix. 

Notice that perturbation theory is valid, as expected, for high energy states determined 
by the expression (3.7). For the high energy spectrum the harmonic oscillator potential acts 
as a small perturbation. Thus the first-order corrections in the energy and the wave function 



are small. Fig. 3.7 shows that the main component of the 105th exact wave function comes 



from the 105th ID box wave function, as it should for small perturbations. 

For low energy states, perturbation theory around the ID box states is not appropriate 
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Figure 3.7: Non-zero components of the 105th exact eigenvector in the basis of a free particle 
in a one-dimensional box. Parameters of the Hamiltonian are u = 16, L = 7r/2, h = m = 1. 
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Figure 3.8: Non-zero components of the third harmonic oscillator eigenvector as expanded 
in the basis of a free particle in a one-dimensional box. Parameters of the Hamiltonian are 
u = IQ, L = 7r/2, h = m = 1. 



since the harmonic oscillator states are the true states in this region. Specifically, for m = 
h = 2L/tt = 1 and u = 16, the first ten states are exactly the harmonic oscillator states 
with a very high accuracy. The next ten states have high overlaps with the corresponding 
harmonic oscillator wave functions. For example, starting from 0.999999 at the tenth state, 
the overlaps go down to 0.880755 at the twentieth state; after that the overlaps get small 
very quickly. Fig. 3.8 shows the structure of the third exact eigenvector when expanded in 
the ID box basis. Notice that the third ID box wave function is almost missing from the 
structure of the third harmonic oscillator wave function. Such small overlap can happen at 
particular values of the parameter uL"^ relevant for the problem at hand. 

This pattern of having a small component of the exact wave function along the corre- 
sponding ID box wave function continues to persist into the transition region. This is an 
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unexpected behavior considering the fact that the first order estimates of the energy levels 
are relatively good. Thus we are confronted with a situation where perturbation theory is 
not appropriate since level spacing is smaller than the magnitude of the "perturbing po- 
tential" but the expectation values of the full Hamiltonian are relatively close to the exact 



eigenvalues , even thought the corresponding ID box wave functions are not at all present 



in the exact wave function as shown in Fig. 3.9 



In conclusion, there is a clear shortcut when using the oblique-basis idea. This allows 
one to use the correct wave functions in the relevant low and high energy regimes relative 
to Ec- There is a clear coherent mixing in the transition region. Such a phenomenon has 
been also observed in the lower p/-shell nuclei 44-48rpj ^^^^ ^^Cr which will be discussed in 
the next chapters. Due to the simplicity of the model, there is a small numerical gain in 
using oblique-basis calculations. However, there could be other cases with a significant gain 
in using the oblique-basis type calculations. Nuclear physics provides one such example as 
demonstrated in our study of ^'^Mg [8]. Another two- mode system of interest is a particle 
confined in two dimensions by an external magnetic field. This system is interesting because 
there is a lifting of the infinite degeneracy of Landau like states due to the confinement 



'A simple explanation of this effect is that the unperturbed energies are such that > 6E^ 
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Figure 3.9: Coherent structure with respect to the non-zero components of the 25th, 27th and 
29th exact eigenvector in the basis of a free particle in a one-dimensional box. Parameters 
of the Hamiltonian are u = 16, L = tt /2, h = m = 1. 



32 



Chapter 4 

Oblique Shell-Model Basics 



Some modem shell-model codes are based on the so-called m-scheme logic, namely, the 
model space is spanned by many-particle configurations (Slater determinants) with good 
third component of the total angular momentum (Mj) [52ll53]. A good total angu- 
lar momentum J, which is a conserved symmetry due to the isotropy of the space, is obtained 
either by angular momentum projection before, after, or as a consequence of the diagonaliza- 
tion of the Hamiltonian. Codes of this type normally achieve a good description of nuclear 
phenomena dominated by the single-particle effects. In these codes the basis consists of sin- 

k 

gle machine words representing the many-particle configurations \ni...nk) = Yl i^t)^" 

s=l 

Unfortunately, an equally good description of collective phenomena within the framework of 
this approach is difficult to obtain due to the computational problems associated with the 
size of the needed model space. On the other hand, the S'?7(3)-based shell-model scheme is 
designed to give a simple interpretation of the collective nuclear phenomena. An ideal sce- 
nario would incorporate both the single-particle degree of freedom and the collective degree 
of freedom, allowing the Hamiltonian of the system to "choose" the admixture that is most 
appropriate. In this chapter, we discuss some of the computational methods and techniques 
used in our calculations. 

4.1 Generalized Eigenvalue Problem 

The usual procedure for solving an eigenvalue problem Hv = Xv is to cast it into a matrix 
equation. In a non-orthogonal basis [54J, this matrix equation includes an overlap matrix 
(0jj = and has the form 

J2iH^JV,-\e.,v,)=0. (4.1) 
j 
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For an orthonormal basis the overlap matrix becomes the identity matrix (0jj — i- Sij), and 
the matrix form of the eigenvalue problem is 



ijVj = Xvi. (4.2) 



When the overlap matrix is positive-definite, the Cholesky algorithm [55j, which de- 
composes into the product of an upper diagonal matrix (U) and its transposed (W^), 



—7- UU , can be used to cast the generalized eigenvalue problem (4.1) back into the 



standard matrix equation (4.2): 



H'v' = Xv', H' = U-^H {U-^f , v' = U^v. (4.3) 
The use of the Cholesky algorithm is essential for identifying the linearly dependent vectors 



within the oblique basis. For large spaces, the effective eigenvalue problem (4.3 ) can be solved 
efficiently by using an appropriately modified Lanczos algorithm which we will discuss in a 
later section. 

For the calculations that will be discussed later, we use two basis sets. The first set 
consists of spherical shell-model states (ssm-states) expressed in spherical single-particle 
coordinates (nlj). The second set has a good SU(3) structure (su3-states) which track 
nuclear deformation this basis set is given in cylindrical single-particle coordinates. By 
construction, both sets have the third projection Mj of the total angular momentum J as a 
good quantum number [22l[T9]. Schematically, these basis vectors and their overlap matrix 
can be represented in the following way: 



basis vectors 



overlap matrix 



Hamiltonian matrix 







ssm 
su3 

1 



basis 
basis 



1 



H 



H, 
H. 



ssmx ssm 



suSxssm 



Hssmx su3 
HguS Xsu3 



Hal3 
Hii3 



(4.4) 
(4.5) 
(4.6) 



In the above, a and i span the following ranges: a = 1,..., dim(ssm-basis) and i = 1,..., 
dim(su3-basis). 

Calculations in a nonorthogonal oblique basis require an evaluation of the matrix elements 
of physical operators plus a knowledge of the scalar product (e^ ■ Ei) related to the overlap 
matrix. While it may be desirable to have an analytical expression for the overlap matrix, as 
we have for the single-particle overlap matrix [57] , for practical purposes it suffices either to 
know the representation of each basis state in a common set that spans the full space, which 
is counter to the overall objective of reducing the number of basis states to a manageable 
subset, or to expand one set in terms of the other. For the present work, the e^, which can 
be represented by a single machine word in a spherical single-particle scheme, were expanded 
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in a cylindrical basis, which is the representation for our collective SU(3) basis vectors. This 
transformation is handled by an efficient routine that exploits two computational aids: bit 
manipulation via logical operations and a weighted search tree for fast data storage and 
retrieval [18]. Transformation of this type has to be done at least once per ssm basis state 
(ca). We transform the ssm-basis states since the result is usually a vector with fewer 
components than a typical SU(3) basis state. There is a simple way to calculate the overlap 
between states in different single-particle bases [58]. However, for the calculation of matrix 
elements of the Hamiltonian, it is better to transform each vector in the basis used by 
the SU(3) states. 

Matrix elements of the one-body and two-body Hamiltonian 



have to be evaluated in each subspace {Hap and Hji), as well as between the two spaces {Ha 



and Hj^), see (4.6). The Ha/s part is normally given and evaluated in a spherical single- 
particle basis. By transforming the Hamiltonian to a cylindrical single-particle basis one can 
obtain the Hji part of H. In order to compute the off-diagonal blocks {Hai, Hjp, and overlap 
matrix elements between SU(3) and ssm-basis states), both basis sets are expanded in a basis 
of Slater determinants using cylindrical single-particle states. For example, any vector within 
the two irreps (8,4) and (9,2) of ^'^Mg has at most 2120 cylindrical Slater determinants; each 
ssm state, which itself is a single spherical Slater determinant, typically expands into less than 
1296 cylindrical Slater determinants. We do not expand the SU(3) states into spherical-basis 
Slater determinants because that would require a significant fraction of the entire spherical 
shell-model space, defeating the rationale of our approach. Taking into account the significant 
number of Hamiltonian matrix elements {Hij and -ffj/j) between multi-component states, it 
should be clear that this is the most time consuming part of the calculation. The extra labels 
associated with the intrinsic quadrupole moment e of each basis state is used to produce well- 
structured band-like matrices and to speed up the calculation. Specifically, basis states are 
pre-ordered according to their deformation as reflected by e, and during the evaluation of H 
a Ae selection rule is applied. 

It is important to point out that knowledge of the overlap matrix O and the matrix 
elements of H in the two spaces {Hajs, Hij) is not enough to obtain the correct off-diagonal 
block Hai- This is clear from the following explicit expression for Hai which contains a 
summation along (/3) that lies outside of the oblique model space 

Hai = ^ HapQpi + ^ HapQpi- 

Thus a direct evaluation of Hai is required. 
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4.2 Geometrical Visualization of the Oblique Basis 



It is instructive to consider a geometrical visualization of the oblique-basis concept. Since 
a set of vectors defines a hyperplane, it is natural to ask the question: "What is the angle 
between hyperplanes defined by the bases under consideration?" To answer this question, 
first consider the angle 9 between a normalized SU(3) basis vector v and the subspace V 
spanned by the spherical shell-model basis vectors. The length of the projected vector vy &V 
is given by cos(w, V) = cos 6* = \vv\- The space V of the spherical shell-model basis vectors 
induces a natural basis in the SU(3) space (n^ = n\Ei). The angle between each new 
basis vector and the space V will again be the length of its projection into the space V, 
but it has the nice property that this set of orthogonal basis vectors stays orthogonal after 
the projection into the space V: 

cos 6 1, = 
riev = 

\nev\ = 

In matrix notation this reads 

where the natural basis vectors are eigenvectors of the symmetric matrix B ■ G"^ 

e . . = e^He. (4.7) 

It follows that IrTevI^ = rie ■ G ■ ■ = e'^n^ ■ = e^, and thus the matrix Q ■ is positive 
definite (|n£y|^ = £^ > 0) with eigenvalues determined by the cos6'. This construction allows 
for a simple visualization of the space spanned by the oblique basis: Choose the x-axis 
to correspond to the space V of all the spherical shell-model basis vectors and represent 
the SU(3) space as a collection of unit vectors each at an angle cos6' = e with respect to 
the X-axis. This construction will be applied later to the geometry of oblique-basis space 
calculations to demonstrate the relative orthogonality of the two vector sets, and Ei. 

4.3 The Lanczos Algorithm 

The Lanczos algorithm is an essential scheme for obtaining a small number of eigenvectors 
corresponding to the lowest or highest eigenvalues [201 EH] • It has been applied successfully 
to spatial dimensions on the order of 10^ and even pushed up to 10^ [59j. This algorithm is 
a simple and very efficient method to build a basis of the Hilbert space associated with an 
eigenvalue problem for an operator H. In its simplest form, one starts with a trial state and 
applies H over and over to generate new states; the process can be applied as many times 
as desired. This way one generates an orthonormal basis in which the corresponding matrix 



cos(fie, V) = \n^v \ , 

^nl{Ei ■ ea)ea = ^n*Gi„e„, 

5^(5^<G..)2 = ^<G,„niG,, 



\nev\ = ■ G ■ G ■ n^. 
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of the operator H is tri-diagonal. The method is recursive and could be used in numerical, 
as well as in analytic calculations [60] . For our toy model described earlier, we have used 
analytic realization (coordinate representation) of the algorithm while for the calculations in 
nuclei a numerical matrix realization has been more suitable due to the Fock representation 
of the states. 

In brief, the algorithm starts with the choice of a first normalized vector vi {{vi\vi) = 1). 
In matrix calculations this vector is often chosen randomly. Then H is applied on vi and a 
new vector orthogonal to Vi is constructed, V2 = Hvi — {vi\Hvi) vi. Next V2 is normalized and 
used to generate a new vector and so on. It can be shown that the basis {vn} generated this 
way is orthonormal, and H is tri-diagonal in this basis. However, numerical noise destroys 
the orthogonality and requires one to do full reorthogonalization of each newly generated 
vector with all previous vectors. One important feature of the Lanczos algorithm is that 
at each new iteration it provides the vector that has the next most important contribution 
within the model space. However, this is true only if our first vector is a good trial guess to 
an exact state. A trial vector with some bad components can cause problems. 

Another feature of the Lanczos algorithm is that it preserves the symmetry of the initial 
vector when this symmetry is a symmetry of the Hamiltonian as wellj^ For example, if H 
is invariant under parity transformation, then the algorithm will produce only even parity 
states. Similarly, if we start with a state with good J and Mj, then applications of H 
will only produce states of the same symmetry. While this can be viewed as an advantage, 
sometimes it can be a problem especially when the symmetry has a finite irreducible sub- 
space. For example, if we start with a vector from a finite irreducible sub-space, then after a 
finite number of iterations the algorithm will exhaust the sub-space and any new vector will 
be linearly dependent on the previously generated vectors. This breakdown of the algorithm 
is generally overcome by introducing a new guess vector. In large matrix diagonalizations, 
the new vector could be a random Gaussian vector, or it could be "the next vector" from 
a prior given set. In our toy model, we carried out Lanczos type calculations using the 
harmonic-oscillator basis as a prior given set. 

An interesting variation to the Lanczos algorithm aiming at a particular k-th eigenvector 
has been suggested by Davidson [61]. In the Davidson algorithm, one tries to increase the 
speed of convergence by modifying the way a new vector is generated. For example, the 
Lanczos algorithm uses Wn = Hvn as a seed for a new vector Vn+i while the Davidson 
algorithm uses the vector Wn = (AfcJ — diag {H))~^ {Hck — Xk^k) with enhanced components 
along the k-th eigenvector. In the Davidson expression for Wn, Ck and are the approximate 
k-th eigenvector and eigenvalue after the n-th iteration, and diag (H) is the diagonal part of 
H. 



To solve the generalized eigenvalue problem (4.3), we have to modify the Lanczos algo- 
rithm. In doing so, it is more efficient to perform consecutive action of the matrices U and H 
on the resulting vectors. The computational time in this case grows like the dimensionality 



^The Lanczos algorithm should in principle conserve symmetries; however, machine round-ofF error often 
mixes in different states. The round-off error is also the reason for performing a complete re-orthogonalization 
for each newly generated state. 
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of the space to the second power (n^). This is to be compared to the case when one would 
first fully multiply these three matrices: (f/~^)^, H, and which grows as the third power 
of n (n^), and then act on vectors. 

In closing this section, we would like to point out some possible future applications of 
the Lanczos algorithm and its modifications. Although recently the algorithm is mostly used 
in huge matrix diagonalizations W^ . it can also have some applications to constrained 
problems as mentioned in our toy model discussion of a particle in a box in the harmonic- 
oscillator basis. For such problems, one would have to project m„ = Hvn in the space 
determined by the constraints, and only after that proceed with the calculation. Another 
interesting application is related to the long standing problem of doing ab-initio calculations 
in nuclear physics with effective interactions derived from a NN-interaction. Some current 
advances in this field which takes advantage of the Lanczos algorithm has been reported by 
Haxton et al. [63| IM]. A link between the Lanczos method and space projection techniques, 
such as Brueckner, Feshbach, and Bloch-Horowitz projection treatments [36l [65l [66] also 
seems intriguing [63] . 

4.4 Mixed- Symmetry Basis for the Nuclear Shell Model 

Even though the purpose of this work has been reiterated several times in different con- 
texts, we feel strong motivation to state it again, but this time in the context of a pure 
practical curiosity. Specifically, can we do calculations using two or more important ba- 
sis sets as usually employed in different shell models? In particular, can we use spherical 
shell-model states, which are related to the single-particle j-shell symmetry (®^-''^^f/(l) = 
U{1) (g) ... (g) U{1) ® f/(l)), together with SU{3) shell-model states, which are related to the 
Q ■ Q interaction? With these questions in mind, which will be answered in the positive sense 
in the next chapter, we continue our examination of the oblique basis concept by focusing 
on the structure of the two basis sets used in our mixed-symmetry shell-model calculations. 
In the rest of the chapter, we briefly go over the spherical single-particle basis, which is then 
followed by a discussion of the SU{3) symmetry-adapted basis. 

4.4.1 Spherical Basis for Single-Particle Excitations 

We have already introduced the main concepts and notations in a previous chapter, as 
well as the basic idea for configuration truncation in the m-scheme basis. Here, we would 
like to touch upon some details related to the specifics of our calculations. First of all the 
spherical basis states are taken from an old m-scheme code NUCK (GLASGOW) [21]. This 
code has been used as a benchmark and a testing ground of our oblique code results. The 
output from NUCK containing the spherical shell-model states (ssm-states) is used as input 
for GlsgwBasis2Redstick which gives a binary form of the proton-neutron basis states for use 
by the oblique code su3pn. In order to speed the calculations when SU{3) states are also 
included, there is an option to use ^-sorted ssm-states (sorting is provided by EpsSorting). 
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4.4.2 SU(3) Basis for Collective Excitations 

We have already reviewed Elliott's SU{3) shell model [4] and some of the single-particle 
labeling schemes in a previous chapter. In this section, we briefly discuss the SU (3) package 
used to generate SU{3) symmetry-adapted states with good third component of the total 
angular momentum. Next the SU{3) single-particle shell-model basis states are discussed, 
and the action of the SU (3) generators is explained. It is then followed by a discussion of the 
structure of the Extreme Weight State(s) (EWS) of an SU{3) irrep for protons (neutrons), 
and especially the Highest Weight State(s) (HWS) of the so-called leading SU{?>) irreps. 
Once a HWS is known, then all states of the corresponding irrep can be constructed using 
SU{3) step operators [23]. Proton (neutron) states with good third component of the angular 
momentum (Mj) are obtained by considering the direct product SU{?>) ® SUs{2) using spin 
highest weight states, as well as spin Lowest Weight States(s) (LWS) [22]. Once the proton 
and neutron highest weight states are constructed, then they can be coupled to different 
possible proton-neutron highest weight states: 

{SU{3) ® SUs{2)), ® {SU{3) ® SUs{2))n ^ {SU{3) ® SUs{2))pn. 

This is the so-called strong coupling scheme. This scheme is used to couple the proton 
and neutron irreps to some final proton-neutron irreps. Each extreme weight state can be 
used for the generation of all states of good Mj within a given SU (3) irrep. Another possible 
coupling scheme extends SUs{2) to S'f/(4) D SUs{2) ® SUt{2). This is the supermultiplet 
scheme which is a good symmetry for light nuclei. Since the goodness of the supermultiplet 
scheme does not extend to heavy nuclei, we will not considered it further in this study. 

• The SU(3) Basis Generator 

The package for generation of SU (3) symmetry-adapted states with good third component 
of the total angular momentum consists of two major codes: (1) SU (3) proton-neutron HWS 
generator {SUS^HWS^GEN) and (2) proton-neutron generator of good Mj (PNGGMJ). The 
SU3.HWS.GEN routine provides the input for PNGGMJ which generates the basis states 
needed for our oblique calculations. 

The overall algorithm has four basic components: 1) definition of the single-particle levels 
and matrix elements of the SU{3) generators for a given proton (neutron) shell; 2) generation 
of the HWS of SU{3) (g) SUs{2) for a given spin 5* and number of protons (neutrons) N; 3) 
coupling of the proton HWS and neutron HWS to the desirable proton-neutron SU{3) HWS 
and SUs{2) LWS; 4) generation of all proton-neutron SU {3)®SUs{2) states with good third 
component of the total angular, Mj. 

The PNGGMJ code accepts any pn-HWS and generates all the states with a given Mj 
value. However, the SUS^HWS^GEN code does not generate all possible pn-HWS since 
one purpose of the current project is to include a few essential SU{3) basis states in an m- 
scheme type calculation. Thus the SU3.HWS.GEN code is set to generate only the leading 
proton and neutron configurations and their possible couplings. Therefore, an additional 
code {SU3Lister) is needed to allow for a quick look at all the proton-neutron SU{3) irreps 
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that can be generated by the current version of the SU3.HWS.GEN code. As a reference 
to a complete hst of proton-neutron SU{?)) irreps, one can use the SU (3) reduced matrix 
elements package [TTj. If other HWS are desired, then the SUS^HWS^GEN code can be 
modified to generate non-leading HWS. This could be done either by using Bahri's method 
[TTj or perhaps by using another more direct approach tailored to the particular application 
at hand. The generation of non-leading HWS is important when one wishes to include states 
that are not maximally deformed in their intrinsic configuration. For example, non-leading 
HWS will be important if the non-Q ■ Q parts of the interaction play a significant role. 



• Single-Particle Levels and the SU(3) Matrix Elements 

The foundation of a microscopic type symmetry-adapted shell-model calculation is the 
structure of the single-particle levels (SPL). The single-particle levels should be related to a 
representation of the symmetry group which is SU{'i) in our case. Therefore, a discussion 
focused on the SU{?)) single-particle levels and matrix elements of the SU (3) generators is 
desirable and will be given in the following few paragraphs. 



Single-particle Levels - Ordering Scheme: Single-particle levels of the r] = 0,1,2,... 

sd, ...) harmonic oscillator shell belong to the symmetric (r^, 0) irrep of S'f/(3). Because 
/i = 0, a typical three-dimensional representation of SU{?>) basis states. Fig. 4.1, reduces to 
a special two-dimensional triangular shape {e and Hp become linearly dependent). Fig 



4.2 



Also, because SU{^) is a compact group, its irreps are finite dimensional, and many-particle 
(fermion) configurations can be conveniently represented as binary strings with a 1 or 
symbolizing the presence or absence of a particle in the corresponding single-particle level. 
(The latter, together with a "sign rule" to accommodate fermion statistics, is a convenient 
computer implementation of a Slater determinant representation of the basis states.) 

Recall that the canonical reduction, S'f/(3) D U{1) ® SU{2), has two additive labels e 
(Qo) and m/ (Lq), and the allowed values of these labels for fixed SU{3) irrep (A,/i) are given 

by m- 



e = 2X + fi - 3{p + q), Up = fi + [p - q), mi = Hp -2m (4.8) 

where 0<p<A, 0<g<;U, and < m < Up. 

A convenient ordering scheme (which tracks the arrows in Fig. 4.2) is set by requiring a 
simple representation of many-particle configurations with maximum quadrupole deforma- 
tion. This objective can be achieved if the states are ordered by e (quadrupole moment) first 
and then by mi (third component of the angular momentum). 



Action of SU(3) Generators on Single-Particle States: To be able to apply the SU{?>) 
generators on many-particle configurations, it suffices to know the action of these generators 
on the single-particle states. The eight generators of 5^7(3) belong to the self-adjoint (1, 1) 
irrep of SU (3). The operator structure should be chosen in the most convenient form for the 
application under consideration. For the present application, this choice is the same as used 
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Figure 4.1: Three-dimensional view of the (A,/i) SU{3) irrep together with the action of the 
SU (3) step operators. 
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sd- shell r\=2 




Figure 4.2: Ordering scheme of the single-particle levels. 



for the basis states, namely, the SU{3) D SU{2) ® U{1) reduction. The matrix elements 
of the S'[/(3) generators can be obtained either by using an application of the appropriate 
Wigner-Eckart theorem or by using explicit expressions [23] for determining the action of 
the operators on the basis states. For computational purposes, it is better to adopt a direct 
solution. We use the fact that the action of the operators is on a product of single-particle 
levels each of which belongs to a symmetric (77,0) irrep of SU{?>). This allows the matrix 
elements of the SU (3) generators to be calculated using properties of the SU (2) only (Fig. 



43j). 

A key feature is the fact that the six non-diagonal generators of S'?7(3) (recall that Lq 
and Qo are diagonal) are rising or lowering generators of SU{2) subgroups of S'[/(3). The 



three SU{2) subgroups and their respective actions are shown in Fig. 4.3 States that are 



collinear with one of the sides of the triangular shape shown on the right in the figure form 
an irrep of the corresponding SU{2) subgroup. 



• Action of SU(3) Generators on Many-Particle States 

Having been supplied with the single-particle levels and the action of the S'?7(3) gener- 
ators on them, we can construct many-particle states and extend the action of the 5'f/(3) 
generators to these many-particle states as well. Since one goal of the oblique-basis project 
is to include essential SU (3) basis states in an m-scheme type calculation, we would briefly 
discuss the maximally deformed HWS for protons (neutrons). These HWS are the leading 
proton (neutron) irreps and can be coupled easily to the leading and other proton-neutron 
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Figure 4.3: Action of the SU{3) generators on the single particle levels. The diagram on 
the left shows that applying an SU (3) generator to a single-particle state results in another 
single-particle state. The diagram on the right shows the action of the six non-diagonal 
generators of SU{3). The vertical solid lines represent the action of the SU{2) subgroup 
that enters the SU{3) D SU{2) ® f/(l) chain. 
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irreps for a given nucleus. Once we have an SU{3) HWS ^SUs{2) LWS state, we can easily 
generate all the states with good Mj within this SU{3) SUs{2) irrep. 



Highest Weight States of SU(3)(8)SU5(2) for Leading Irreps: So far we have con- 
structed single-particle states and evaluated matrix elements of the generators of SU{3) when 
they act on these states. The next step is to construct many-particle HWS of SU {3)®SUsi2). 
In the chosen SU{3) labelling scheme, there are seven extreme states which correspond to 



the vertex points of the three-dimensional diagram (Fig. 4.1) of a general (A,/i) irrep. We 



are particularly interested in the vertex that has the maximum value for the quadrupole 



moment of the system (Fig. 4.1). Our HWS is the state with e = 2X + /j,, Up = n, and 
mi = fi. This HWS (maximum value of m/ for maximum e) can be easily constructed by 
ensuring that the action of the SU{3) rising generators annihilates it. Indeed, for such a 
HWS, the values of A and can be determined from its e and mi labels. 

Selecting the leading (A,/i) irrep (HWS with maximum overall value of e) out of all 
possible irreps of an N fermion system with total system spin S is very simple within the 
chosen scheme. This is because the number of particles with spin up and spin down is 
uniquely determined by the solution of two linear equations: 

N = n^ + 72^, 2S = n^- n^. 

The second of these two equations expresses the fact that we also require the state to be 
highest weight with respect to SUs{2). Further, maximizing the value of Qo is achieved by 



filling the single-particle states of the {rj^O) irrep (Fig. 4.2) from bottom to top. The chosen 
scheme ensures that this simple procedure gives maximum values for e and mi. The SU{3) 
irrep labels (A, /i) are obtained by evaluating the quadrupole moment (Qo) and angular 
momentum projection (Lq) which are additive quantum numbers. 

For example, in the S(i-shell, there are six single-particle levels corresponding to the (2, 0) 
irrep of SU{3). The HWS of the leading irrep for = 6 particles and total system spin 



^ = 1 is (3,3) (Fig. 4.4), whereas for 5 = the leading SU{3) irrep is (6,0) (Fig. 4.5). 



These many-particle configurations are HWS with respect to SU{3) and SUs{2). 

Generating SU(3) States with Step Operators: Once we have the HWS of SU{3) 
® SUs{2), we can generate any other state of the SU{3) irrep (A,/i) by applying step 
operators similar to those given by Hecht. This process is needed when we produce proton- 



neutron coupled SU{3) irreps. By using the parameterization (4.8) [23], we can identify the 
corresponding step operators. Fig. |4.1 



It is important to note that applying p-move or g-move step operators to states on the top 
surface yields other states (or zero) on that same surface. Since the states on the top surface 
are HWS with respect to SU{2) in the SU{3) D SU{2) ® U{1) reduction, the m-move step 
operator is an SU{2) lowering operator which changes the third component of the angular 
momentum (mi). The p-move and g-move step operators can be obtained by imposing the 
restriction that they generate only transformations within the SU{2) HWS space. From an 
algebraic perspective, the p-move and the m-move operators are linear in SU{3) generators 
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Figure 4.4: Highest weight state of the leading (3, 3) irrep for = 6 and 5 = 1 in the 
S(i-shell. 
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Figure 4.5: Highest weight state of the leading (6, 0) irrep for = 6 and S* = in the 
sd-sheW. 
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while the g-move operator is quadratic. Nevertheless, the state generation process can be 
written in such a way that the g-move operator effectively reduces to a linear action. These 
step operators can also be obtained by a projection operator technique |67j . 



Generating Proton-Neutron SU(3) Highest Weight States: In brief, the generation 
of the proton-neutron SU{3) HWS is just a matter of SU (3) and SUs{2) couplings. However, 
the actual algorithm for the generation process is somewhat backwards to the structure of the 
sentences that we would use to describe it. What we mean by this is that when the number 
of protons and the number of neutrons are given together with their harmonic oscillator 
shells, then the algorithm does a loop over all possible total proton-neutron spins (Spn). For 
each total spin (Spn), loops are made over the possible proton spin (Sp) and neutron spin 
(Sn) that can couple to the total spin (S'p„). This way the proton HWS can be constructed, 
as described in the previous sections, by using the proton number(A'p) and spin (Sp) to get 
the {SU{3) ® SUs{2))p HWS state. The same is done for the neutron HWS. This way all 
the major labels, \N{X, fj,)e,np,mi)p ® \SMs)p for protons and \N{X, ^)e,np,mi)p ® \SMs)p 
for neutrons, have been determined and one can proceed with the details of the coupling to 

{SU{3) ® SUs{2))p ® iSU{3) ® SUs{2))^ ^ {SU{3) ® SUs{2))p,,. 

One final detail on the generation process is that the proton HWS state is actually an 
SU{3) HWS but a spin LWS while for the neutrons it is an 5*^7(3) and a spin HWS. Then 
they are coupled to a proton-neutron SU (3) HWS with a spin LWS structure. The reason 
behind this is that in the SU{2) type couphng one does a loop such that nip + m„ = mp„. 
Hence, if irip is the minimal m- value of the proton spin such that irtp + m„ = irtpn for a fixed 
mpni then the m„ should be the maximal m-value of the neutron spin. Therefore, the loop 
which satisfies nip + rUn = nipn will have an increase of rup and simultaneous decrease of m„ 
so that rUpn stays fixed. This is also the reason why the total proton-neutron state is a LWS 
{Ms = —Spn) SO that the coupling to the good Mj = Mi + Ms is done in the same way since 
Ml is at maximum in the proton-neutron SU{3) HWS. 



Generating States of Good Mj: Since the action of the SU{3) commutes with that of 
the spin group {SUs{2)), it is not difficult to achieve the final goal of states with good third 
component of the total angular momentum (Mj). Recall that we have just generated proton- 
neutron EWS {SU{3)^^^ (8) SUs{2)l\ys)- We also introduced the procedure to generate 
other states of an SU{3) irrep by applying step operators. Each of these states remains a 
LWS of SUs{2). Hence, after each move on the top surface (see Fig. 4.1), using p-move 
or g-move step operators, we can apply spin rising and angular momentum lowering to the 
corresponding state toward Mj = M / + Ms. This way, we can generate all proton-neutron 
states with labels: N, S, (A,/i), e, Up, Mi, and Mj. 
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Chapter 5 

^^Mg Mixed- Symmetry Calculations 



The success and applicability of the oblique-basis approach to ^^Mg, which will be demon- 
strated in the following sections, can be related to the fact that the spherical shell-model 
states are eigenstates of the one-body Hamiltonian (^£ja^aj) while the two-body part of 
the Hamiltonian j Vki,ijO'to,'j'akai) is strongly correlated with the quadrupole-quadrupole 
interaction {Q ■ Q) which is diagonal in the SU(3) basis [68j- By combining spherical shell- 
model states and SU(3) states, one accommodates, from the onset, the dominant modes of 
the system. 

In this chapter we discuss the oblique-basis technique as applied to ^^Mg [8]. This is 
a strongly deformed nucleus with well-known collective properties and is one of the best 
manifestations of the Elliott's SU(3) symmetry [4]. In terms of dimensionality of the model 
space, adding a few leading SU(3) irreps to a highly truncated spherical shell-model basis 
results in significant gains in the convergence of the low-energy spectra towards the full space 
result. In particular, the addition of leading SU(3) irreps yields the right placement of the 
K=2 band and the correct order for most of the low-lying levels. Indeed, an even more 
detailed analysis shows that the structure of the low-lying states is significantly improved 
through the addition of a few SU(3) irreps. The Hamiltonian used in our analysis is the 
Wildenthal interaction [12j. 

In the following sections we summarize some of the important features of the spherical, 
SU(3), and mixed-symmetry type shell-model calculations. First we discuss the dimensions 
of each model space. Then we consider the ground-state energy as a function of the model 
space used. Next we focus on the structure of the low-energy spectrum, and finally we discuss 
the structure of the states as compared to the exact sd-shell results. 

5.1 Structure of the Model Spaces 

One important question in any computational study is the dimensions of the matrices 
involved, as well as the structure of the model space used. In this section, we address this 
question by briefly summarizing the space structure and dimensions for the spherical, SU(3), 
and oblique shell-model calculations. 
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Table 5.1: Labels and Mj = dimensions for various ^^Mg calculations. The leading SU(3) 
irrep is denoted by (8,4) while (8,4)&(9,2) implies that (9,2) irreps have also been included. 
The SM(n) spaces correspond to spherical shell-model partitions with n valence particles 
excited out of the ^5/2 shell into the si/2 and ^3/2 levels. 

Model space (8,4) (8,4)fc(9,2) SM{0) SM{1) SM{2) gM(4) FULL 
space dimension 23 128 29 449 2829 18290 28503 

% of the full space 0.08 0.45 0.10 1.57 9.92 64.17 100 



5.1.1 Model Space Dimensions 

Our model space for ^"^Mg consists of 4 valence protons and 4 valence neutrons in the 
Qhjjj Sfi-shell. The m-scheme dimensionality (Mj = 0) of this space is 28503. This space is 
denoted as FULL in the figures that follow. To test the effects of truncations, calculations 
were also carried out permitting n particles to be excited out of the lowest orbit, i.e. 
d^2^d2,/2Si/2Y, and are denoted as SM(n). The SM(2) approximation is of particular interest 
since it allows one to take into account the effect of pairing correlations (one pair maximum) 
in the 'secondary levels' {si/2 and 0^3/2 for the ^"^Mg) with a minimum expansion of the 
model space. The SU(3) part of the basis includes two scenarios: one with only the leading 
representation of SU(3), which for ^^Mg is the (8,4) irrep, with dimensionality 23 for the 
Mj = space and denoted in what follows by (8,4); and another with the (8,4) irrep plus 
the next most important representation of SU(3), namely the (9,2). The (9,2) irrep occurs 
three times, once with S = {Mj = dimensionality 15) and twice with S = 1 {Mj = 
dimensionality 2 x 45 = 90). All three (9,2) irreps have total Mj = dimensionality of 
15+90=105. The (8,4)&(9,2) case has total Mj = dimensionality of 23+105=128 and is 



denoted by (8,4)&(9,2). In Table 5.1 we summarize the dimensionalities involved. 



5.1.2 Visualizing the Oblique Basis for ^^Mg 

After obtaining an idea of the space dimensions involved, we now try to visualize the 
oblique basis. The method described in the previous chapter can be used to visualize the 
structure of the oblique basis. First, consider the SM(2) space enhanced by the SU(3) irreps 
(8,4)&(9,2). Since the SM(2) and (8,4)&(9,2) spaces are both relatively small (see Table 



5.1), we expect the basis vectors of these spaces to be nearly orthogonal. This orthogonality 



is clearly seen from inset (a) in Fig. |5.1 Inset (b) in Fig. 5.1 shows a loss of orthogonality 
between the SM(4) and the (8,4)&(9,2) basis vectors. This is due to the fact that SM(4) 
space is about 64% of the full sd-space. Therefore, there is a relatively high probability 
that some linear combinations of SU(3) basis vectors lie in the SM(4) space. Indeed, it 
can be shown that there are five vectors from (8,4)&;(9,2) that lie within the SM(4) space. 
Such redundant vectors are identified and excluded from the calculation within the Cholesky 
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SU(3) basis V^^" SU(3) basis 




SM(2) space 




Figure 5.1: Orthogonality of the basis vectors in the obhque geometry. The SU(3) space 
consists of (8,4)&(9,2) basis vectors with the shell-model spaces (SM(n) with n=2 and 4) 
indicated by a horizontal line, (a) SM(2) and the natural SU(3) basis vectors and (b) SM(4) 
and the natural SU(3) basis vectors. In the latter case (b), there are five SU(3) vectors that 
lie in the SM(4) space. 

algorithm when it is applied to the overlap matrix (O — ?■ UU^). 

5.2 Spectral Characteristics 

Reproducing the correct energy spectra of a nucleus is one of the goals of any nuclear- 
structure study. Since we are trying to develop a new concept for nuclear structure stud- 
ies, the mixed-symmetry approach, we are currently comparing our results only with full 
shell-model calculations. Therefore, a computation-to-computation comparison is our real- 
ity check. In the next sub-sections we compare the ground-state energy and energy spec- 
trum for ^^Mg as calculated with the Wildenthal interaction [12] using spherical, SU(3), and 
mixed-symmetry shell-model bases. 

5.2.1 Ground-State Energy 

We now turn to the consideration of the main results of the oblique-basis calculation. 



starting with ground-state convergence issues. The results shown in Fig. |5.2| illustrate 
that the oblique-basis calculation gives good dimensional convergence in the sense that the 
calculated ground-state energy for the SM(2)-|-(8,4)&(9,2) calculation is 3.3 MeV below the 
calculated energy for the SM(2) space alone. Adding the SU(3) irreps only increases the 
size of the space from 9.9% to 10.4% of the full space. Compare this 0.5% increase in the 
size of the space with the huge (54%) increase in going from SM(2) to an SM(4) calculation. 
For the latter, the ground-state energy is 4.2 MeV lower than the SM(2) result, somewhat 
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better than for the SM(2) + (8,4)&(9,2) calculation but in 64.2% rather than 10.4% of the 
full model space. 

The exponential fall-off of the ground-state energy in Fig. 



5.2 is striking. It has been 



observed many times, and has been recently suggested as a possible extrapolation procedure 
for obtaining the ground-state energy [02] • An even more rigorous extrapolation procedure 
has been suggested by Mizusaki and Imada [TD]. Within this procedure, one can also estimate 
the error of a given calculation. Their procedure is also applicable to other observables, as 
well. 



5.2.2 Energy and Angular Momentum of the Low Lying States 



Fig. 5.3 and Fig. 5.4 show that the oblique-basis calculation positions the K=2 band head 



correctly. Furthermore, most of the other low-energy levels are also positioned correctly. The 



results for pure spherical and pure SU(3) calculations are shown in Fig. 5.3[ As can be seen 
from the results in Fig. 

notice that in 



5.3, an SM(4) calculation (64% of the full model space) is needed 

Also, 



to get the ordering of the lowest angular momentum states correctly, 
this case the third and fourth energy levels are practically degenerate. On the other hand, 
it only takes 0.5% of the full space to achieve comparable success with SU(3). In particular. 
Fig. 5.3 shows that an SU(3) calculation using only the (8,4) and (9,2) irreps gives the 



right ordering of the lowest levels. Note that the first few low-energy levels for SM(2) are 
close in energy to the corresponding low-energy levels for the (8,4)&;(9,2) result. Since these 



two spaces are nearly orthogonal (see Fig. 5.1 ), these two sets of levels mix strongly in an 



oblique calculation and yield excellent results. The comparable ground-state energies of the 



SM(2) and (8,4)&(9,2) configurations can also be seen in Fig. 5.2 



Compare the spectra shown in Fig. |5.3| with the results from the oblique-basis calculations 
shown in Fig. 5.4 From this comparison one can see that the correct level structure can 
be achieved by using 1.6% (SM(l)-|-(8,4)) of the full sd-space. However, one should also 
notice that for the SM(0) + (8,4) space, which is only 0.2% of the full space and the minimum 
obhque-basis calculation, the results are quite close to the correct level structure. Despite 
the fact that the ground-state energy of the oblique-basis calculation is higher than the 
ground-state energy for the SM(4)-type calculation, the oblique calculations are favorable in 
terms of dimensionality considerations and correctness of the level structure. 



5.3 Overlaps with the Full 5(i-shell Calculation 



Figs. |5.5| - |5.8| focus on the actual structure of the states by showing overlaps of eigenstates 
calculated in the SM(n), SU(3), and oblique bases with the corresponding states of the full 



space calculation. Specifically, in Fig. 5.5, overlaps of states for pure SM(n) and pure SU(3)- 



type calculations are given. Note that the SM(4) states have big overlap (90%) for the first 
few eigenstates. This should not be too surprising since SM(4) covers 64% of the full space. 

The results in Fig. 5.5 show that in general SU(3)-based calculations give much better 
results than the low- dimensional SM(n)-type calculations. The SM(n)-based calculations 
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Figure 5.2: Calculated ground-state energy for ^''Mg. Ground-state energy as a function 
of the various model spaces. Note the dramatic increase in binding (3.3 MeV) in going 
from SM(2) to SM(2)+two SU(3) irreps, (8,4)&(9,2), (a 0.5% increase in the dimensionality 
of the model space). Enlarging the space from SM(2) to SM(4) (a 54% increase in the 
dimensionality of the model space) adds 4.2 MeV in binding energy. 
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Figure 5.3: Structure of the energy levels for Mg for different calculations. Pure m-scheme 
spherical basis calculations are on the left-hand side of the graph; pure SU(3) basis calcu- 
lations are on the right-hand side; the spectrum from the FULL space calculation is in the 
center. 
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Figure 5.4: Energy levels for ^'^Mg as calculated for different oblique bases. The SM(4) basis 
calculation is included for comparison. 
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Figure 5.5: Overlaps of the pure spherical and SU(3) with the FULL states. The first four 
bars represent the SM(0), SM(1), SM(2), and SM(4) calculations, the next two bars represent 
SU(3) calculations, etc. 
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have irregular overlaps along the low- lying states and require SM(4), which is 64% of the 
full space, to get relatively well behaved overlaps. This can be seen most clearly from the 
inset labeled SM in Fig. 5.6 Note that the SM(0) contributions to the third, fifth, and sixth 
states are very low while SM(1) and SM(2) have varying contributions. The structure of the 
SU(3)-type states leads to a stable picture for the oblique calculations as shown in the inset 
SM(n) + (8,4) and SM(n) + (8,4)&(9,2) in Fig. ^ 
In Fig. 

From this graph, one can see that the 



5.7[ the improvement in the structure of the calculated states is followed as 

basis. 



the SU(3) states are added to the SM(n 
improvement to the SM(0)- and SM(l)-type calculation is due mainly to the goodness of 
SU(3) itself. The improvement obtained in the oblique calculation is due to the SU(3) 
enhancement of the SM(2) space. From this graph, one can also conclude that there is only 
a small gain in going to the SM(4)-based oblique calculation. However, this improvement 
can not be achieved by any other means with such a small increase in the model space. This 
is clear from a careful examination of Fig. 5.2 where one can see that the SM(5) result. 



which has 25142 basis vectors (88% of the full sci-space), gives the same ground-state energy 
as the SM(4) + (8,4)&(9,2) resuh (64.6% of the full sd-space). 

Finally, to compare the three schemes - SU(3), SM(n), and the various oblique-basis 
combinations - representative overlaps are shown in Fig. |5.8[ From these results, it is very 
clear that SU(3)-type basis states yield the right structure in a very low order. In particular, 
in Fig. 5.8, it can be seen that a 90% overlap with the exact eigenvectors can be achieved 
by using only 10% of the total space, SM(2) + (8,4)&(9,2). Furthermore, Fig. 5.8 also shows 
that SU(3) enhances the SM(4) results yielding eigenstates with overlaps that are very close 
(~ 98%) to the exact results. 
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Figure 5.6: Overlaps of the oblique-basis states with the exact states (set 1). Inset SM 
contains the overlaps for the pure spherical shell-model basis states only. Inset SM-|-(8,4) 
contains the overlaps of the SM basis enhanced by the leading SU(3) irrep (8,4). Inset 
SM+(8,4)&(9,2) has the (9,2) irreps included as well. 
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Figure 5.7: Overlaps of the oblique-basis states with the exact states (set II). Each inset 
represents a particular SM(n)-type calculation, showing how the overlaps change along the 
corresponding oblique-basis calculation. 
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Figure 5.8: Representative overlaps of pure SM(n), pure SU(3), and oblique-basis results 
with the exact full srf-shell eigenstates. A number within a bar denotes the state with the 
overlap shown by the bar if it is different from the number for the exact full-space calculation 
shown on the abscissa. For example, for SM(2) the third eigenvector overlaps the most with 
the fourth exact eigenstate, not the third, while the fifth SM(2) eigenvector has the overlap 
shown with the third exact eigenstate. 
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Chapter 6 

Study of Lower p/-shell Nuclei 



For ^^Mg, the single-particle excitations, described by the spherical shell model, and the 
collective excitations, described by the SU{3) shell model, are of comparable importance. 
In the previous chapter, we have shown the relevance of the oblique-basis calculation for 
^^Mg. It is, therefore, natural to seek other nuclear systems to apply the mixed-symmetry 
method to. The even-even nuclei in the sd-shell are one place to start. For the sd-shell 
nuclei, however, one can perform full S(i-shell calculations with modern computer codes. 

The p/-shell nuclei are another option. For these nuclei, full p/-shell calculations have 
just recently been achieved [59l|72]. However, adding only the leading and next to the leading 
irreps, as it has been done for ^^Mg, is not sufficient for the lower p/-shell nuclei, Ti and 
Cr, to obtain results as good as those for ^^Mg. This is because the spherical shell model 
provides a significant part of the low-energy wave functions of these nuclei within a few 
spherical shell-model configurations, while in the SU{3) shell-model basis one needs more 
than a few SU{3) irreps. This is due mainly to the strong breaking of SU (3) in the lower 
p/-shell induced by the spin-orbit interaction fl]. When the spin-orbit splitting is removed, 
the importance of the SU (3) basis is restored. Although the usual SU (3) structure of the 
states is lost, there is an adiabatic SU{3) mixing which gives rise to the coherent structure of 
the yrast states. We have already seen this coherent mixing phenomenon in our toy model. 
In nuclei, however, this coherent mixing can be interpreted as an illustration of the intrinsic 
state idea where all the states within a given band can be projected out from a particular 
intrinsic state Even more, in nuclei this coherent structure is assumed to be a result 
of an adiabatic SU (3) mixing which means that some observables stay close to the SU (3) 
limit, that is, as if there was a pure SU{3) symmetry. Specifically, the B{E2) values remain 
strongly enhanced with values close to the SU{3) symmetry limit. It is important to point 
out that there is a coherent mixing of the spherical shell-model states as well. 

In this chapter, we will discuss our study of the even-even lower p/-shell nuclei 44-48rjj 
and ^'^Cr. First, we show what a few SU{3) irreps can do for us within a mixed-symmetry 
calculations for ^^Ti nucleus. Then, we focus on the spin-orbit interaction that strongly 
breaks the SU{3) symmetry. We conclude this chapter by discussing the coherent structure 
of the states and the adiabatic SU{3) mixing which produces enhanced B{E2) values. 
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Table 6.1: Labels and Mj = dimensions for various "^^Ti oblique calculations. 



Model space (12,0) (12, 0)fc(10, 1) SM{0) SM{1) SM{2) SM{A) FULL 
space dimension 7 84 72 580 1908 3360 4000 

% of the full space 0.18 2.1 1.8 14.5 47.7 84 100 



6.1 Ti Oblique-Basis Results 

The simplest even-even nucleus in the p/-shell, from a computational point of view, is 
^^Ti. In this section we discuss our oblique-basis calculations for ^'^Ti. If one compares 
the spherical shell-model with the SU{?)) shell-model results within the framework of a 
realistic interaction, such as KB3 interaction |9], then S'f/(3) seems to be badly broken. 
Specifically, the ground-state energy and wave function are poorly reproduced. This seems 
to be a common trend in the even-even srf-shell nuclei as well [71]. Even the ground-state 
energies within the oblique-basis calculations do not look prominent. However, at a closer 
examination one finds that the oblique-basis idea still works. The results may be not as good 
as in ^^Mg, but there are some close analogies. For example, the SM(1) space in ^^Ti seems 
to be what SM(2) is for ^'^Mg, while the SM(2) space in ^^Ti seems to be what SM(4) is for 
^■^Mg. By that we mean that the SU{?)) enhanced SM(1) basis in "^^Ti gives overlaps that 
are compatible with the overlaps of the pure SM(2) calculation. In the next few sub-sections 
we briefly illustrate these findings. 



6.1.1 Model Space Dimensions 

The model space for '^^Ti consists of 2 valence protons and 2 valence neutrons in the pf- 



shell. We use the same notation for the m-scheme spherical bases as in 5.1 The SU(3) 
part of the basis includes the leading irrep of SU(3), which for ^^Ti is (12,0) with Mj = 
dimensionality 7, and the next to the leading irrep, namely the (10,1). The (10,1) irrep occurs 
three times, once with S = {Mj = dimensionality 11) and twice with S = 1 {Mj = 
dimensionality 2 x 33 = 66). All three (10,1) irreps have total Mj = dimensionality of 
11+66=77. The (12,0)&(10,1) case has total Mj = dimensionality of 7+77=84 and is 



denoted by (12,0)&(10,1). In Table 6.1 we summarize the dimensionalities involved. 

As in the case of ^^Mg, there are linearly dependent vectors within some of the ^'^Ti 
calculations. For example, there is one such vector in the SM(2) + (12,0) space, two in the 
SM(3)+(12,0), two in the SM(1) + (12,0)&(10,1), twelve in the SM(2)-F(12,0)&(10,1), and 
thirty-three in the SM(3) + (12,0)&(10,1). Each linearly dependent vector is handled as 
discussed in the chapter devoted to ^^Mg. 
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Figure 6.1: Ground-state energy for ^^Ti as a function of the various model spaces. The 
SU(3) irreps used are (12,0) and (10,1). 



6.1.2 Ground-State Energy 

Fig. 6.1 shows that the obhque-basis calculation of the ground-state energy for ^'^Ti 
does not give results as good as for ^^Mg. The calculated ground-state energy for the 
SM(1)+(12,0)&(10,1) calculation is 0.85 MeV below the calculated energy for the SM(1) 
space alone. Adding the two SU(3) irreps to the SM(1) space increases the size of the space 
from 14.5% to 16.6% of the full space. This is a 2.1% increase, while going from the SM(1) 
to SM(2) involves an increase of 33.2% in the model space. For the latter, the ground-state 
energy is 2.2 MeV lower than the SM(1) result. 



6.1.3 Energy Spectrum of the Low Lying States 

We have seen that the position of the K=2 band head for ^^Mg is correct for the SU(3)- 
type calculations but not for the low-dimensional SM(n) calculations. In ^"^Ti this seems to 
be the opposite, the SM(n)-type calculations reproduce the position of the K=2 band head 



while the SU(3) do not, as shown in the upper graph in Fig. 6.2 Furthermore, most of the 
low-energy levels are much higher for the pure SU(3) limit than for the pure SM(n) case. 



61 



Thus, one may expect that these two sets of levels (the SM(n) and the SU(3)) may not mix 
as strongly in an oblique calculation as for the ^^Mg case. Surprisingly, the oblique-basis 
calculations seem to produce a good spectral structure as shown in the lower graph in Fig. 



6.2, Notice that the SM(2)-|-(12,0)&(10,1) spectrum is very good and compatible with the 



SM(3) spectrum. This is 50% compared to 84% of the model space. 



6.1.4 Overlaps with the Exact Calculation 



The top graph in Fig. 6.3 shows overlaps of states for pure SM(n) and pure SU(3)-type 
calculations while the lower part shows some selected overlaps from the oblique calculations. 
Notice that the overlaps of the pure SU(3)-type calculations are very small, often less than 
40%, while the SM(n) results are far better with the SM(2)-type calculations having about 
80% overlap with the exact states. Note that the SM(3) states have big overlap (>97%) for 
the first few eigenstates. This is not surprising since SM(3) covers 84% of the full space. 
What is surprisingly good is that the SM(2)-|-(12,0)&(10,l)-type calculation is as good as 
the SM(3). Even more, the SM(1) + (12,0)&(10,1) overlaps seem to be often bigger than the 
SM(2) overlaps. 



6.2 Set Up for the Study of the SU{3) Breaking 

To understand better the results of the mixed-mode calculations described in the previous 
section, we need to recall that the oblique-basis method is expected to work well if we are 
dealing with two or more competing and compatible modes. Therefore, if the Hamiltonian 
of the system is dominated by its one-body term, then the effect of the two-body part will be 
suppressed. However, if the single-particle energies are degenerate, the importance of SU{3) 
should reappear. In the next subsections, we discuss the structure of the Hamiltonian, as 
well as some of the computational methods used in our calculations. 



6.2.1 Interaction Hamiltonian 

To retain clarity of the discussion, we recall the structure and notations of the one- plus 
two-body Hamiltonian: 

H = ^Eia^ai + ^ Vki^ija^a'^akai. 

The summation indexes range over the single-particle levels included in the model space. 
We only consider levels of the p/-shell which have the following radial (n), orbital (/) and 
total angular momentum (j) quantum numbers: nlj = 10/7/2,0/5/2,1^3/2,1^1/2}- In what 
follows, the radial quantum number (n) is dropped since the Ij labels provide a unique 
labelling scheme for single-shell applications. It is common practice to replace the four 
single-particle energies Ei by the P and I ■ s interactions: 'Ylii^i'^t'^i ^ ^i^i " ^i^i ' ~ 
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Figure 6.2: Structure of the energy levels for ^^Ti for different calculations. Pure m-scheme 
spherical-basis calculations are on the left-hand side of the upper graph; pure SU(3)-basis 
calculations are on the right-hand side; the spectrum from the FULL space calculation is in 
the center. The spectra form oblique-basis calculations are in the lower graph. 
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Figure 6.3: ^^Ti wave function overlaps of pure spherical, SU(3), and oblique states with the 
FULL states. The first four bars in the upper graph represent the SM(0), SM(1), SM(2), and 
SM(3) calculations, the next two bars represent the SU(3) calculations, etc. Representative 
overlaps of pure SM(n), pure SU(3), and oblique-basis results with the exact full p/-shell 
eigenstates are shown in the lower graph. 
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where e is the average binding energy per valence particle, rii counts the total number of 
valence particles, and a and P are dimensionless parameters giving the interaction strength 
of the P and / ■ s terms. For realistic single-particle energies used in the KB3 interaction 



(6.1), these parameters are e = 2.6 MeV, (5 = 0.0096, = 1.3333, and a/ = 1.7143. The 
small value of (5 signals small splitting (6.2 ) while the values of a demonstrate the presence 



of a strong spin-orbit splitting. 

A significant part of the two-body interaction, Vki^ij, maps onto the quadrupole-quadrupole 
{Q ■ Q) and the pairing (P) interactions. Since Q ■ Q can be written in terms of SU{?>) gener- 
ators, it induces no SU{i) breaking, as has been discussed in the first chapters. Hence Q ■ Q 
serves to re-enforce the importance of the Elliott model [1] . However, the pairing interaction 
mixes different SU{?>) irreps, but in our study it does not seem to cause any strong SU{?>) 
breaking. In this analysis the two-body part of the Hamiltonian {Vki^ij ) is fixed by the Kuo- 
Brown-3 (KB3) interaction matrix elements, and the single-particle energies, e^, are changed 
as described below. 

The following single-particle energies are normally used with the KB3 interaction [9]: 

KB3 [MeV] : Sp^ = 4, Sp, =2, Sf = 6, = 0. (6.1) 

3 3 f 2 

For the purposes of the current study, it is important to know the centroids of the p- and 
/-shells. For example, the energy centroid of the p-shell is given by: 

dim(pi) + dim(p3) 

7 2 2 

Sr. 



'P 



dim(pi) + dim(p3) 



In what follows, we label by KBSp.f that Hamiltonian which uses the KB3 two-body 
interaction with single-particle p- and /-shell energies set to their centroid values: 

KB3p_/ [MeV] : Sp, = Sp^ = 2.6670, €f = Sj, = 2.5710. (6.2) 

3 3 3 2 

We use KBSpf for the case when the single-particle energies are set to their overall 
average: 



KB3p/ [MeV] : £p = = 2.6 (6.3) 
Due to the near degeneracy of the single-particle energies of the KB3p-f interaction 



(6.2), the results for the KB3pf case are very similar to those for KB?>p_f. 



6.2.2 Computational Procedures 

In our study, we have focused on ^^Ti, ^^Ti, "^^Ti, and ^'^Cr because these are p/-shell 
equivalents of ^°Ne, ^^Ne, ^^Ne, and ^^Mg, respectively, which are known to be good SU{?)) 
sd-shell nuclei. Furthermore, data on these nuclei are readily available from the National 
Nuclear Data Center (NNDC) [73] and full p/-shell calculations are feasible [72] . The model 
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Table 6.2: Space dimensions for the m-scheme calculations in the full p/-shell model space. 
We have used even parity and even isospin basis states with no restrictions on the total 
angular momentum J except for the Mj = case where only states with even J values have 
been selected. 

Nucleus Mj = Mr = 6 Mr = 10 Mr = U 



134 
3846 
8997 



Ti 1080 514 30 

'Ti 43630 32297 4693 

'Ti 317972 278610 57876 

Cr 492724 451857 104658 



dimensionalities for full-space calculations increase very rapidly when approaching the mid- 



shell region; those for the cases considered here are given in Table 6.2 



The computational procedures and tools used in the analysis of the SU(3) symmetry 
breaking are described in this section. In brief, the Hamiltonian and other matrices are 
calculated using an m-scheme shell-model code [12] while the eigenvectors and eigenvalues 
are obtained by means of the Lanczos algorithm [21]. AH the calculations are done in the 
full p/-shell model space. 

First, the Hamiltonian H for each interaction {KB3 (6.1), KBSp^f (6.2), and KB3pf 



(6.3)) is generated. Then the eigenvalues and eigenvectors are calculated and the yrast 
states identified. Next, the matrix for the second-order Casimir operator of SU{3) , namely 
C2 = (3L^ + Q ■Q)/4:, is generated using the shell- model code, and a moments method [H] is 
used to diagonalize the C2 matrix by starting the Lanczos procedure with specific eigenvectors 
of H for which an 5*^7(3) decomposition is desired. Finally, B{E2) values in e^fm^ units 
are calculated from one-body densities using Siegert's theorem with a typical value for the 
effective charge [SHI ES]' ?e// = 0.5, so = (1 -|- q'e//)e = 1.5e and e„ = (q'e//)e = 0.5e. 

Even though the procedure can generate the spectral decomposition of a state in terms 
of the eigenvectors of C2 of S't/(3), this alone is not sufficient to determine uniquely all irrep 
labels A and n of SU{?>). For example, C2 has the same eigenvalue for the (A,/x) and A) 
irreps. Nevertheless, since for the first few leading irreps (largest C2 values) the A and /x 
values can be uniquely determined [76], this procedure suffices for our study. 

Usually, when considering full-space calculations, a balance between computer time and 
accuracy has to be considered. While the Lanczos algorithm [21] is known to yield a good 
approximation for the lowest or highest eigenvalues and eigenvectors, it normally does a 
relatively poor job for intermediate states. This means, for example, that higher states, in 
particular high total angular momentum states, may be poorly represented or, in a worst case 
scenario, not show up at all when these states are close to or beyond the truncation edge of 
the chosen submatrix. An obvious way to maintain a good approximation is to run the code 
for each Mj value, that is, Mj = 0, 2, 4, 6. . . . However, this might be a very time consuming 
process, but nonetheless one which could be reduced significantly if only a few Mj values 
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are used for each run. For the calculations of this study, we have used Mj = 0, 6, 10, and 
14. To maintain high confidence in the approximation of the intermediate states which have 
J = 2,4,8,12,... we required that they be within the first half of all the states produced. 
The code output was set for 29 states. A further verification of the accuracy of the procedure 
is whether the energies of the same state calculated using different Mj runs are close to one 
another. For example, as a consistency check the energy of the lowest J = 6 state in the 
Mj = run was compared to the energy of the same state obtained from the Mj = 6 run. 



6.3 Measuring Symmetry Breaking Using C2 of SU{3) 

In this section, we discuss the results of our study on the SU (3) symmetry breaking in the 
p/-shell. In order to identify the SU{3) structure of an yrast state, we calculate the spectral 
distribution of the state along the second Casimir operator (C2) of SU{3) as described in 
the previous section. From the spectral distribution, we can clearly determine whether the 
5*^7(3) symmetry is broken or not. However, a graphic or table representation of the data 
becomes very inelegant with growing space dimensions. Thus, we have decided to use also 
average quantities, such as centroid, width, and skewness of the distributions to illustrate 
the main points one can deduce from a complicated spectral distribution. 



6.3.1 Spectral Distribution 



The first set, Figs-6.4 and 6.5, demonstrates the recovery of the SU{3) symmetry as the 
single-particle spin-orbit interaction is turned off, that is, in going from the KB3 to the 
KB3p-f interaction. Corresponding results for the KB3pf interaction are similar to the 
KB3p-f results. In each graph, C2 values of SU{3) are given on the horizontal axis with 
the contribution of each SU (3) state on the vertical axis. The bars within each cluster are 
contributions to the yrast states starting with the ground state (J = 0) on the left. Hence 
the second bar in each cluster is for the J = 2 yrast state, etc. 

We have chosen ''^Ti for an in-depth consideration of the fragmentation of the C2 strength 



in yrast states. The results for the nondegenerate KB3 interaction are shown in Fig. 6.4 
In this case the highest contribution (biggest bar) is more than 50% which corresponds to a 
C2 value of 114 for the J = 12 state. The C2 = 114 value is for (A,/i) = (8,2) which is two 
SU{3) irreps down from the leading one, (A,/i) = (12,0) with C2 = 180. The leading irrep 
only contributes about 10% to the J = 12 yrast state. The contribution of the next to the 
leading irrep, C2 = 144 for (A, /x) = (10, 1), is slightly less than 40%. Thus, for all practical 
purposes, the first three irreps determine the structure of the J = 12 yrast state. This 
illustrates that the high total angular momentum J states are composed of only the first few 
SU (3) irreps. This is easily understood because high J values require high orbital angular 
momentum (L) states which are only present in SU{3) irreps with large C2 values. The high 
J states may therefore be considered to be states with good SU (3) symmetry. However, this 
is not the case with the ground state of ^^Ti which has very important contributions from 
states with C2 values 60, 72, 90, 114, 144, and 180 with respective percentages, 7.5, 25, 10, 
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Figure 6.4: Strength distribution of C2 of SU(3) in yrast states of ^^Ti for realistic single- 
particle energies with Kuo-Brown-3 two-body interaction {KB3). 
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Figure 6.5: Strength distribution of C2 of SU (3) in yrast states of '^^Ti for degenerate single- 
particle energies with Kuo-Brown-3 two-body interaction (KBSp^f). 



69 



21, 8, and 21%. This shows that the leading irrep is not the biggest contributor to the J = 
ground state; there are two other contributors with about 20%, the third {C2 = 114) and 
seventh (C2 = 72) SU{3) irrep. 

When the spin-orbit interaction is turned off, which yields nearly degenerate single- 
particle energies since the single-particle orbit-orbit splitting is small, one has the KBSp.f 
interaction, and in this case the structure of the yrast states changes dramatically, as shown 



in Fig. 6^ In Fig. 6^ one can see that the leading irrep plays a dominant role as its 
contribution is now more than 50% of every yrast state. As in the previous case, the high 
total angular momentum J states have the biggest contributions from the leading irrep, for 
example, more than 97% for J = 12, 91% for J = 10, and 80% for J = 8. The ground state 
is composed of few irreps with C2 values 72, 114, and 180, but in this case the leading irrep 
with C2 = 180 makes up more than 52% of the total with the other two most important 
irreps contributing 21% [C2 = 72, (A,/i) = (4,4)] and 23% [C2 = 114, (A,/i) = (8,2)]. 

6.3.2 Moments of the Spectral Distributions. 



An alternative way to show the recovery of the SU{3) symmetry is given in Fig. 6.6 and 



Fig. |6.7[ These figures show the centroid, width, and skewness of the C2 distributions. 
The J values are plotted on the horizontal axis with the centroids given on the vertical 
axis. The width of the distribution is indicated by the length of the error bars which is 

just the rms deviation, AC2 = \J {{C2 — (6*2))^), from the average value of the second-order 



Casimir operator (C2). The third central moment, bC2 = ^ {{C2 - (^2))^, which measures 

the asymmetry, is indicated by the length of the error bar above, AC2 + and below, 
AC2 — the average value. 

Note that the recovery of the leading irrep when the spin-orbit interaction is turned off is 
clearly signaled not only through an increase in the absolute values of the first centroid (C2) 
but also through the skewness 5C2- For example, in ^"^Ti with the KB?) interaction (spin- 
orbit interaction turned on) the ground state J = has (C2) = 110 and skewness 5C2 = 33. 
This changes for the KBSp.f interaction to (C2) = 139 and a skewness of 6C2 = —37, as 



shown in Fig. 6.6 The equivalent of the ^^Ti graph for the ^^Ti case is shown in Fig. 6.7 
As for the ^^Ti case, the results show the recovery of the SU (3) symmetry in ^^Ti when the 
single-particle spin-orbit interaction is turned off. 

6.3.3 Coherent Spectral Structure 

We now turn to a discussion of the coherence nature of the yrast states. First, notice 



that the widths of the distributions as defined by AC2 = y {{C2 — {C2)) ) are surprisingly 



unaffected (Fig. 6.6 and Fig. 6.7) when turning the spin-orbit interaction on and off. This 



effect occurs in all cases studied: ^^Ti, ^^Ti, ^^Ti, and ^^Cr. The more detailed graphs. Fig. 



6.4 and Fig. 6.5, offer an explanation in terms of the fragmentation of the C2 distribution. 



As can be seen from these graphs, the irreps that are present in the structure of a given yrast 
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Figure 6.6: Average C2 values for KB3 and KBSp.f interactions in ^*Ti 
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Figure 6.7: Average C2 values for KlSS and KBSp f interactions in ^^Ti. 



state in the presence of the spin-orbit interaction (Fig. 6.4) remain present, even though 
with red uced strength, in the structure of the state when the spin-orbit interaction is turned 

off (Fig. 6.5). As a consequence, AC2 = \J{{C2 - {C2)f) which measures the overall spread 
of contributing irreps, is more or less independent of the spin-orbit interaction. One can see 
a sharp decrease in the width of the distribution only for high spin states like J = 12 in the 
graph for ^^Ti in Fig. 6.6 



Fig. 6.8 demonstrates the coherent nature of the states within the yrast band. The three 



graphs shown give the spectrum of the second-order Casimir operator C2 of SU (3) for the 
J = 0, 2 and 4 yrast states in ^^Cr. The axes are labelled the same way as in Figs. |6.4 and 



6.5 , but in this case all bars are for a single yrast state. In this figure there are three peaks 
surrounded by smaller bars that yield a very similar enveloping shape for the given yrast 
states. The fragmentation and spread of C2 values is nearly identical for these states with 
no dominant irrep, indicative of severe S'f/(3) symmetry breaking. 

Graphs for the KBSp^f case, when the spin-orbit interaction is turned off, are not shown 
since the results are similar to the results for '^^Ti shown in Fig. 6.5 For example, when 



the spin-orbit interaction is on (KB3), the leading irrep for '^^Cr has a C2 value of 396 and 



this accounts for only around 10% of the total strength distribution (see Fig. 6.8), but when 



the spin-orbit interaction is off (KB3p_/), the leading irrep is the dominant irrep with more 
than 55% of the total strength. 

We conclude the section with a discussion of the coherent structure of the yrast states 
by an illustration of the coherent structure of the ^^Cr states within the spherical shell- 
model basis. The inset (a) of Fig. 6.9 shows the spectral structure of the lowest yrast 
states (J = 0,2,4, and 6), as calculated with the KB3 interaction, with respect to the 
spherical configuration basis. Notice the common spectral distribution of these states. The 
distribution along the energy configurations, related to excitation energies smaller than the 
harmonic oscillator spacing {< Ihw = 10 MeV), provides an illustration and support of the 
energy-based configuration truncation scheme. The bump at 12 MeV is probably related to 
the fact that this is only a p/-shell calculation (Ohw) which does not include the multi-shell 
excitations that are at energies above Ihu = 10 MeV. 



6.3.4 Enhanced Electromagnetic Transitions. 

Our results on the lower p/-shell nuclei, so far, have shown that SU{3) symmetry breaking 
in this region is driven by the single-particle spin-orbit splitting. However, even though 
states of the yrast band exhibit SU{3) symmetry breaking, the yrast band B{E2) values 
are insensitive to this fragmentation of the SU{3) symmetry; specifically, the quadrupole 
collectivity as measured by B{E2) transition strengths between low- lying members of the 
yrast band remain high even though SU{3) appears to be broken. 



Relative B{E2) values are shown in Figs. 6.10, 6.11, and 6.12, that is, B{E2) strengths 
normalized to the B{E2 : 2"*" — )■ O"*") value. For isoscalar transitions, the relative B{E2) 
strengths are insensitive to the chosen effective charges which may be used to bring the the- 
oretical B{E2 : 2+ — 0^) numbers into agreement with the experimental values. Whenever 
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Figure 6.8: Coherent structure of the first three yrast states in Cr calculated using realistic 
single-particle energies with Kuo-Brown-3 two-body interaction {KB3). On the horizontal 
axis is C2 of SU (3) with contribution of each SU (3) state to the corresponding yrast state 
on the vertical axis. 
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Figure 6.9: Coherent mixing and SU(3) breaking and recovery in ^^Cr. Inset (a) demon- 
strates the coherent structure of the yrast states with respect to the spherical shell-model 
configuration basis (KB3); (b) coherent structure of the yrast states with respect to the 
SU(3) basis (KB3); (c) recovery of the SU(3) symmetry within the KBSpf interaction. 
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Figure 6.10: Relative B{E2) values ( 1(^2+^0+) ) ^^^i. The 5(E2 : 2+ ^ 0+) transi- 
tion values are 122.69e^/m^ for the experiment, 104.82e^/m^ for the KB3 interaction, and 
138.58e2/rn^ for the KBSp.f case. 



absolute B{E2 : 2+ — )■ 0^) values are given, they are in fm^ units and the effective charges 
are 1.5e for protons and 0.5e for neutrons (q'e// = 0.5). 

The first graph on relative B{E2) values (Fig. 6.10 ) recaps our results for '^^Ti. Cal- 
culated relative B{E2) values for ^^Ti corresponding to the spin-orbit interaction turned on 
(KB3) and spin-orbit interaction off (KB3p_/) are very close to the pure SU{3) limit. The 
agreement with experiment is very satisfactory except for the 4+ — t- 2+ and 8^ — t- 6+ tran- 



sitions. However, the experimental data [73] on the 8^ — t- 6+ transition give only an upper 
limit of 0.5 pico-seconds to the half-life. We have used the worse case, namely a half-life of 
0.5 ps, as a smaller value would increase the relative B{E2). For example, a half-life of 0.05 
ps will agree well with the relative B{E2) value for the KB3p-f interaction. This example 
supports the adiabatic mixing which seems to be present for all the yrast states of ^'^Ti. 

Fig. 6.11 shows B{E2) values for ^^Ti. In this case there are deviations from adiabatic 
mixing for the 6"*" — )■ 4"'', 10^ 8^, and higher transitions. Two experimental data sets are 
shown in Fig. 6.11 data from the NNDC is denoted as Exp_(NNDC), and updated data 
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Figure 6.11: Relative B{E2) values ( 1^2:2+^0+) ) ^^^i. The B{E2 : 2+ ^ 0+) transition 
values are 199.82e^/m^ for the experimental data, ISLTOe^/m"^ for the updated experimental 
data, 208eV^^ for KB3 interaction, and 299.83e^fm'^ for KB3p_f. 



on 2+ — )■ 0^ and 4+ — )■ 2+ transitions from [77\ is denoted as Exp_(Updated). For ^^Ti 
the agreement with the experiment is not as good as for ^^Ti. However, the experimental 
situation is also less certain. The adiabatic behavior is well demonstrated for the first three 
yrast states 0^, 2+, and 4"*" via relative B{E2) values for the KB3 and KBSp.f interactions 
which are very close to the SU (3) limit. 

We conclude this section by showing the recovery of the SU (3) symmetry; this time via 
relative B{E2) values as shown for ^'^Ti in Fig. 6.12, In Fig. 6.12 we see that for the 



degenerate single particle case (KB3p_/) the first few transitions have relative B{E2) values 
which follow the SU{3) limit very closely. On the other hand, the interaction involving 
spin-orbit splitting (KB3) is far from the 5*^7(3) limit. The B{E2 : 4+ — )■ 2"*") transition is 
strongly enhanced due to the adiabatic mixing which is missing in the higher than J = 4 
yrast states. 
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Figure 6.12: Relative B(E2) values (^^ir^^) for ^^Ti. The B(E2 : 2+ ^ 0+) transition 
values are 144.23e^/m^ for the experimental data, 155.5e^/m^ for the updated experimental 
data, 202.4eV^^ for KB3 interaction, and 445.32eV^^ for KBSpJ. 
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Chapter 7 



Summary and Discussions 



The primary goal of the current work has been to study and apply a new method-the 
mixed-symmetry approach-for large shell-model calculations. Our aim was to combine two 
very successful computational methods: the m-scheme spherical shell model and the 5*^7(3) 
shell model. In the process of this study, we have realized a new computational paradigm: an 
oblique-basis calculation that can be used to capture the mixed-mode structure of complex 
systems, such as the atomic nuclei. 

The two methods, the m-scheme and SU{3), are closely connected to the two dominant 
but often competing modes that characterize the structure of atomic nuclei: the single- 
particle shell structure underpinned by the validity of the mean-field concept, and the many- 
particle collective behavior manifested through the nuclear deformation. This is reflected in 
two dominant elements in the nuclear Hamiltonian: the single-particle term. Hi = Y^-EiUi, 
and a collective two-body term i^2- The collective term H2 is dominated by the quadrupole- 
quadrupole interaction, Hqq = Q ■ Q which has good SU{3) symmetry. It follows that the 
simplified Hamiltonian H = Yli ^i^i ~ xQ ' Q has two exactly solvable limits and thus can 
be considered to be a two-mode system. 

To probe the nature of such a system, we have considered a simple toy model: the one- 
dimensional harmonic oscillator in a box. As for real nuclei, this system has a finite volume 
and a restoring force whose potential is of a harmonic oscillator type. For this model, there 
is a well-defined energy scale which measures the strength of the potential at the boundary 
of the box, Ec = For this system, the use of two sets of basis vectors, one for each of 

the two limits, has physical appeal, especially at energies near Ef.. One basis set consists of 
the harmonic oscillator states; the other set consists of basis states of a particle in a box. In 
the regime of strong mixing of the two modes at an energy scale compatible with Ec, there 
is a coherent structure expressed through a quasi-perturbational behavior of the system. 
Specifically, in this energy region first-order perturbation theory is not appropriate since the 
zeroth order approximation to the wave function is very poor; nevertheless, the first-order 
estimates of the energies are very close to the actual results. Even more, the structure of 
the exact wave functions exhibits a coherent mixing (Fig. 3.9) similar to the one observed 
in nuclei (Fig. 6.8). 
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An application of the mixed-symmetry basis calculations to ^^Mg, using the realistic 
USD interaction of Wildenthal, has served to demonstrate the validity of the mixed-mode 
shell-model scheme. In this case, the oblique-basis consists of the traditional spherical states, 
which yield a diagonal representation of the single-particle interaction, together with collec- 
tive SU(3) configurations, which yield a diagonal quadrupole-quadrupole interaction. The 
results obtained in a space that spans less than 10% of the full-space reproduce the correct 
binding energy, within 2% of the full-space result, as well as the low-energy spectrum, and 
the structure of the states within 90% overlap with the exact states. In contrast, for an 
m-scheme spherical shell-model calculation, one needs about 60% of the full space to obtain 
results comparable with the oblique basis results. Calculations for '^''Ti also support the 
mixed-mode shell-model scheme, even though calculations using a few SU{3) irreps are not 
as good as the standard spherical shell-model calculations. And, as the results confirmed, 
the combined basis yields less enhancements. For example, an oblique-basis calculation in 
50% of the full p/-shell space is as good as a usual m-scheme calculation in 80% of space. 
These results show very clearly that if the important modes can be isolated, then one can 
build an oblique theory that incorporates leading configurations of each mode and could get 
good convergence in a limited model space. 

The study of the lower p/-shell nuclei 44-48'jj ^^^^ '^^Cr, using the realistic Kuo-Brown-3 
(KB3) interaction, has shown strong SU(3) symmetry breaking due mainly to the single- 
particle spin-orbit splitting. When the spin-orbit splitting is reduced, the importance of the 
SU (3) as seen through a growth in the dominance of the leading irrep is restored. Thus 
the KB3 Hamiltonian is at least a two-mode system. This is further supported by the 
behavior of the yrast band B(E2) values that seem to be insensitive to the fragmentation 
of the SU(3) symmetry. Specifically, the quadrupole collectivity as measured by the B(E2) 
strengths remains high even though the SU(3) symmetry is rather badly broken. This has 
been attributed to a quasi-SU(3) symmetry where the observables behave like a pure SU(3) 
symmetry while the true eigenvectors exhibit a strong coherent structure with respect to 
each of the two bases. This has been observed in all yrast states for the "^^Ti case; while 
for the other nuclei studied, this coherence breaks down after the first few yrast states. 
In particular, even though the yrast states are not dominated by a single SU{3) irrep, the 
B{E2 : 4+ — i- 2+) values remain strongly enhanced with values close (usually within 10-20%) 
to the SU{3) symmetry limit. 

From a technical point of view, there are some other possible basis sets to be studied]^ 
For example, one can try to use deformed Nilsson basis states, or a basis set generated from 
a Hartree-Fock type procedure [71] • One can even try simple cylindrical basis states with an 
appropriate procedure to maintain a complete set for good spin quantum numbers. If good 

^In our study, SU(3) is shown to be good due to the 3D harmonic osciUator and the dominance of the Q.Q 
interaction in nuclei. The cyhndrical basis is just the easiest way to construct the SU(3) states and seems 
to be most economical in terms of components. From computational point of view, a good total angular 
momentum (J) and its third component (Mj) for the SU(3) states are essential. However, if one can find 
any other basis set, besides the SU(3)-based one, with good J and Mj, then things may be as good, or even 
better. 
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rotational symmetry is to be sacrificed, then one can try a Lanczos algorithm which keeps 
only the big components during the iteration process. 

Another further development of the theory and its application is a study of other S(i-shell 
nuclei as well as p/-shell nuclei. Such studies will further test the theory and the codes that 
have been developed. In spite of the results in the lower p/-shell, it is expected that in the 
mid-shell region some sort of SU (3) collective structure is important. Thus, the oblique- 
basis calculation may be an important alternative for calculating structure of nuclei, such 
as ^^Fe and ^^Ni. Another possibility is to integrate the oblique basis concept into no-core 
calculations of the type developed by [78]. Such an extension would involve the symplectic 
group for multi-shell correlations rather than just SU(3) [70] . 

An extension of the theory to a multi-mode oblique shell-model calculation is also a 
possibility. An immediate extension of the current scheme might use the eigenvectors of the 
pairing interaction [80] within the Sp(4) algebraic approach to the nuclear structure [HT] . 
together with the collective SU(3) states and spherical shell-model states. Even the three 
exact limits of the IBM [2] can be considered to comprise a three-mode system. 

Further, an even broader extension of the theory would involve a general procedure for 
the identification of dominant modes from any one- and two-body Hamiltonian along with a 
complementary partitioning of the model space into physically relevant subspaces with small 
overlaps. One can then start with eigenstates for an arbitrary subspace and constructively 
improve the results by including corrections from the remaining subspaces. It should be 
possible to do this by keeping only a small set of the calculated lowest energy states at each 
iteration. Hamiltonian-driven basis sets can also be considered. In particular, the method 
may use eigenstates of the very-near closed shell nuclei obtained from a full shell-model 
calculation to form Hamiltonian driven J-pair states for mid-shell nuclei [3]. This type of 
extension would mimic the Interacting Boson Model (IBM) [33] and the so-called broken-pair 
theory [3] . Nonetheless, the real benefit of this approach is expected when the system is far 
away from any exactly solvable limit of the Hamiltonian and the spaces encountered are too 
large to allow for exact calculations. 

In summary, we have studied a new computational method, the oblique-basis method. 
The concept has been applied to a toy model, as well as to some realistic nuclear systems. 
For realistic nuclei, we used spherical and cylindrical single particle states to perform our 
mixed-symmetry calculations. We have studied ^'^Mg in the sii-shell and ^^Ti in the p/-shell 
in a mixed-symmetry basis. For ^^Mg, we have seen very promising results with respect to the 
energy spectra and the structures of the wave functions. When these results are translated 
into model space dimensions, we see that an oblique-basis calculation in 10% of the full 
sd-shell space is as good as a usual m-scheme calculation in 60% of the full srf-shell space. 
For ^^Ti, the results are less pronounced due to the dominance of the one-body over the two- 
body part of the Hamiltonian. However, in model space dimensions, the results for ^^Ti state 
that an oblique-basis calculation in 50% of the full p/-shell space is as good as a usual m- 
scheme calculation in 80% of the full p/-shell space. Through a detailed study of ^^Ti, "^^Ti, 
^^Ti, and '^^Cr in the full p/-shell, we have confirmed the effect of the one-body part of the 
Hamiltonian, that is, the strong SU (3) symmetry breaking is due to the spin-orbit interaction 
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which sphts the single-particle energies. For degenerate single-particle energies, we have seen 
that one recovers the dominance of the leading SU (3) irrep which is consistent with two- 
body interaction dominated by the quadrupole-quadrupole interaction. Along our study, 
we have seen some interesting coherent structures, such as coherent mixing of basis states, 
quasi-perturbative behavior in the toy model, and an enhanced B{E2) strength toward the 
SU{3) limit in nuclei. 

In concision, the main positive outcome of this work is a prove-of-principle of the mixed- 
mode concept. We have shown that such calculations are doable and may yield better results 
and lead to a clearer understanding of complex systems. Problems yet to be solved are re- 
lated mainly to the software package and its development. First of all, a routine for the 
complete generation of SU(3) shell model basis is needed. Basis sets other than the spherical 
shell-model and SU(3) shell-model basis sets are also desirable; some possible basis sets have 
been discussed. Another important software component is a set of commonly used physical 
observables and their matrix elements. The most important improvement, however, is to 
implement error estimate of the final results and possible extrapolation procedure for esti- 
mating the exact energy eigenvalues. Immediate further work should include a concentrated 
study of other sd-shell nuclei, p/-shell nuclei, and multi-shell calculations. Applications to 
atomic and molecular physics are also possible. 
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Appendix A 



On the Wave Function Spread and 
Localization 

In nuclear physics, we often use the three-dimensional harmonic-oscillator (3D HO) poten- 
tial as a zeroth order approximation to the nuclear mean-field potential. This is usually done 
in the center-of-mass coordinate system, assuming that all the nucleons experience the same 

attractive potential, Hq = + mfl'^x'^^ j2]. If we assume the same localization for the 

nucleons, there seems to be a localization paradox since we are dealing with fermions that 
must obey the Pauli exclusion principle. However, this apparent paradox is resolved by us- 
ing many-particle Slater determinant wavef unctions, constructed by filling the single-particle 
levels of the three-dimensional harmonic-oscillator potential. The Slater determinant form 
satisfies the Pauli principle requirements and yields different localization structure for each 
nucleon. 

A harmonic-oscillator potential is appropriate near stable equilibrium where the inter- 
action potential should have a local minimum. If rotational invariance applies, then the 
potential near stable equilibrium should actually be a three-dimensional harmonic-oscillator 
potential. However, because of the Pauli principle, only the closed-shell nuclei have a spher- 
ical shape, other nuclei have non-spherical ground state distributions that can be character- 
ized as oblate, prolate or tri-axial. As a consequence, for non-closed shell nuclei, a deformed 
three-dimensional harmonic-oscillator potential is a more appropriate "mean field". This 
idea is incorporated in the deformed Nilsson model [6]: 

H = h mVl?x^ + emVl^xi + vul + visl ■ s. 

2m 

Here, e is a measure of the deformation when P and / • s provide for the correct shell closures 
and magic numbers in nuclei. 

Large scale numerical calculations usually use basis functions of the three-dimensional 
harmonic oscillator (3D HO). This way the wave function parameter u will not match the 
corresponding parameter Q in the Hamiltonian. For example, = + ^ may be very 
different from cj = Q in super-deformed nuclei. It is therefore interesting to look at the 
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behavior of the fixed-basis calculations with respect to localization and energy scale for the 
one-dimensional harmonic oscillator (ID HO). 

Here we fix the parameters of the ID HO Hamiltonian tohem = Q = h= l. Thus, 
its spectrum is simple: En = n + ^. The basis consists of displaced and scaled harmonic- 
oscillator wave functions, "^niil + 0/^)'^ the ^ = states are squeezed/stretched states, and 
(7 = 1 states are coherent states [52] ■ All the calculations are done using the default settings 
for Mathematica 4.1 [83]. 

Fig. A.l shows the number of the fixed-basis states needed to achieve convergence to the 
10^^ in the ground-state eigenvalue and eigenvector as a function of the parameter u. The 
convergence criteria for the eigenvalues requires two successive eigenvalues to be less than the 
accuracy limit (10~^) apart. The convergence criteria for the eigenvectors is WH'^ — E^^W < 
the accuracy limit. As expected, the eigenvalues converge much earlier than the eigenvectors. 

Fig. A. 2 (a) shows the first few basis states (a; = 0.2), the harmonic-oscillator potential 



(fi = 1), and the true ground-state wave function. Fig. A. 2 (b) shows the calculated ground- 
state wave functions at different dimensions of the basis states with u = 0.2. From these 
graphs, it is clear that when u < Q = 1, Fig. A. 2 (a) and (b), one uses more and more basis 
states to produce the correct wave function behavior within the classically forbidden region. 
When u > fl = 1, Fig. A. 2 (c) and (d), the focus is actually concentrated on getting the 



correct shape of the wave functions within the potential well. 

Fig. |A.3 is similar to Fig. A.l but shows the convergence within the displaced (coherent 
states) harmonic-oscillator wave function basis {"^niQ + O)- -Due to parity conservation, there 
is a good symmetry under the ^ — )■ — ^ transformation. An example of the basis structure 
and convergence path similar to Fig. A.2| is shown in Fig. A. 4 for C, = —2. 

Although one should be able to solve any problem in any arbitrarily chosen orthonormal 
basis, the considerations presented here point to the need for properly modified basis states 
to reduce the model-space dimension and thus to avoid problems due to numerical noise. In 
the process of optimizing the basis set for a particular Hamiltonian, the orthogonality of the 
basis would inevitably be destroyed. In the toy model of a two-mode system, a few possible 
types of basis-state refinements were considered. Originally, the oblique-basis method was 
concerned with two or more basis sets as described in the toy model and in nuclear physics 
applications. However, the idea of the basis refinement can be extended in quite a general 
way as described in the next section. 
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Figure A.l: Ground-state convergence for the harmonic-oscillator problem with Q = 1 using 
squeezed basis states. The number of basis states needed for 10"'' convergence accuracy of the 
ground-state eigenvalue is shown in green squares. The red circles are for the ground-state 
eigenf unction. 
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Figure A. 2: Role of the wave function spread, (a) Stretched basis states with uj = 0.2 within 
the harmonic-oscillator potential ^2 = 1. (b) Consecutive approximations of the harmonic- 
oscillator ground state (red) using the stretched basis states, (c) Squeezed basis states with 
u = 4 within the harmonic-oscillator potential = 1. (d) Consecutive approximations of 
the harmonic-oscillator ground state (red) using the squeezed basis states. 
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Figure A. 3: Ground-state convergence for the harmonic-oscillator problem with $7 = 1 using 
coherent basis states. The number of basis states needed for 10""^ convergence accuracy of 
the ground-state eigenvalue is shown with the green squares. The red circles are for the 
ground-state eigenfunction. 
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Figure A. 4: Role of the wave function localization, (a) Coherent basis states with displace- 
ment C, = —2 within the harmonic-oscillator potential = 1. (b) Consecutive approxima- 
tions of the harmonic-oscillator ground state (red) using the coherent basis states. 
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Appendix B 

Variationally-Improved Basis Method 



The usual fixed-basis method can be derived from the Rayleigh-Ritz variation principle. If 
one considers minimization of E (c) with respect to c for a Hamihonian (H) using the basis 
states (pn {x; u): 

E [c] = \H\ ^) - A(c ■ c - 1), ^{x) = J2 Cn0n (x; UJ) 

n 

of a Hermitian operator Y {u) such that 

Y (u) (pn (x; Uj) = YnCpn [x] , 

then 5E[c\/5c*^ = is equivalent to solving the matrix eigenvalue problem Ylm HnmCm = Ac„, 
where Hnm = (0n l-f^l 0m) , and thus the set of As provides information for the eigenvalues 
of H. If the set {0n} is taken to be non-orthogonal, then we have a generalized eigenvalue 

problem Y,m i^nm - ^fJ'nm) Cm = wherC /inm = (0n|0m)- 

Notice that there is a freedom that we have not yet specified: it is the choice of u and 
Y (u). Here, u is the set of parameters characterizing the Hermitian operator Y (u) . Usually 
one fixes u from experience or by simply applying Rayleigh-Ritz variation with respect to u 
in E [u] = (00 i^) \H\ 00 i^))- Thus, we have an orthonormal basis 0„ with the same u for 
any n. One can try other procedures for fixing u and Y {u) as well [27] . 

All this seems fine as long as the spectrum of H is expected to be similar to the spectrum 
of Y (u) , but what if the potential for Y {u) does not match the "landscape" of the potential 
for H. For example, would the harmonic-oscillator potential wave functions be appropriate 
for solving an anharmonic potential problem or a double-well potential? In principle, one 
should be able to use any basis, but it may be at the expense of long and tedious calculations. 
Therefore, we may try to let a; be a free parameter for different basis functions 0„. Then we 
can find Un by variation of E [cu] = {(j)n (oj) \H\(j)n (c^))- 

Often u is related to the relevant energy scale. If we start with the correct wave function, 
but with the 'wrong' parameters, then clearly a variational approach on the parameters will 
give us the right answer immediately. 
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In the case of the harmonic-oscillator wave functions, one can argue that a calculation 
with a varying u is equivalent to a multi-shell calculation with a fixed u parameter. Since the 
harmonic-oscillator basis is a complete basis, then each function 0„ {x,Un) can be expanded 
in the basis associated with Q, for example: 

(f)n {X, Un) = ^ C^0fc (x, fi) . 
k 

Therefore, 0„ [x, Un) can be viewed as the result of a multi-shell calculation with the harmonic- 
oscillator parameter Q. 

Next, we discuss how in general one can refine any initial basis set so that each basis vector 
in the new and improved basis is the optimal one with respect to the Hamiltonian under 
consideration. Then, instead of refining the basis vectors, one can effectively renormalize the 
parameters in the Hamiltonian. 

The main idea is to optimize each trial vector by applying the Rayleigh-Ritz variational 
principle on [\E', ^4]: 

E[^,A] = {A'^\H\A^), 6E['^,A] = 0, 

where A represents the affine group A in R". An element a of ^ has a rotational component 
r and a translational component t, so that {ax)j = r-^Xj + tj. 
If G is the symmetry group of H, then: 

g-'Hg = H, geG. 

Therefore, the Rayleigh-Ritz variational principle should be applied with respect to the ho- 
mogeneous space M = A/G that excludes the symmetry transformations G. For example, 
a translational symmetry of H means that scaling and rotation are the relevant transforma- 
tions. Since the physical systems usually have rotational and translational symmetry, then 
only scaling is left as a relevant operation for constructing a variationally-improved basis: 

^ (Xi, X„) ^ {SXI, SXn) . 

The transformation |\E') — )■ |v4\I') can be defined to maintain the normalization of the states: 

^' (x) (x) = \/det (r)^ (rx + t) . 

However, this transformation is not a unitary transformation in general, and therefore, it will 
not map orthonormal states into a new set of orthogonal states. Using scaling as a variational 
parameter has been done previously. Specifically, in the context of the confined systems it 
was used by Martin and Cruz to study hydrogen and helium enclosed in a spherical shell 

li J dy = J det (r) dx is used when y = rx + t, then the variationally-improved basis can 
be treated as a renormalization problem for the Hamiltonian H: 

E [^, A] = (A^ \H\ A^) = j det (r) ^* (rx + t)H (p, x) * (rx + t)dx = 
= j ^* (y) H {rp, r-' {y - t)) ^ {y) dy. 
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Rescaling in the above way seems to be related more to the scahng methods in condensed 
matter physics. For example, consider the one-dimensional harmonic oscillator: 

2 2^ 

Then, the equation for the scale parameter s from E s] when ^ (x) — )■ ^/s'^ (sx) is: 

Here, (Q^) = |g2| ^) and Aq = y/{Q^). It is assumed that * is such that (Q) = (P) = 0, 
which means that the localization of the wave function has been selected. Evaluating E [i>,s] 
at the extremum = Aq/Ap gives: 

E [*] = AqAp. 

Finally, using [p,q] = —i ^ AqAp > |, wc find that the minimum of the energy is exactly 
the zero point energy for the harmonic oscillator Eq = ^. Notice that quantum mechanics was 
only used to provide us with a constraint on the fluctuations of the observables q and p; other 
than that, we can consider the system as purely statistical. Thus, different AqAp will give 
different value of E [\E']. Turning this argument around, we would expect AqAp > {n+ 
when Aq and Ap are evaluated in the space of wave functions with n-nodcs. 

Another interesting way to obtain the same result is to use H expressed in terms of the 
operators a"^ and a. Then, by using coherent states as trail wave functions a\z) = z\z) we 
have E [z\ — \z\^ + \ and thus Eq — |. 
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Appendix C 



On the Loss of Hermiticity 



When the choice of the basis is not carried out with appropriate attention, an operator, 
supposedly Hermitian, may acquire a non-hermitian matrix reahzation within this basis. 
For example, a wrong basis may produce a non- Hermitian matrix for the Hamiltonian under 
consideration. Although this is unlikely to be encountered within the finite shell-model 
calculations using an occupation number representation, it is an obstacle when one wishes 
to use a hard core potential and a harmonic-oscillator basis [2]. 

Here we discuss the problem of a free particle in a one-dimensional box in the harmonic- 
oscillator basis. In order to proceed, we notice that the Hilbert space for the harmonic 
oscillator is not quite the same as for the free particle in a one-dimensional box. This is clear 
from the domains of the wave functions. The harmonic-oscillator wave functions are defined 
on the whole real axis R^, when the wave functions for a free particle in a box are defined 
on a finite interval [— This discrepancy is easily fixed by projecting the harmonic- 
oscillator wave functions onto the interval [— L,L], which changes the inner product for the 
wave functions: 



However, in this basis the matrix corresponding to H will be nonhermitian in general. 

To understand the loss of hermiticity, we look at the off-diagonal matrix elements of the 
momentum operator (P = —ih-^): 




f* (x) g (x) dx — )■ 




/* (x) g (x) dx. 



L 



L 




L 



L 




-L 



-L 



L 




-L 
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It is clear from the above expression that the hermiticity will be maintained only when 
all of the basis functions are zer^ at the boundary of the interval [-L,L]. This condition is 
essential for solving exactly the quantization of a free particle in a one-dimensional box. 



^Wave functions with the same value at ± L is a necessary condition; wave functions should be zero at 
± L only for an infinite potential. 
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Appendix D 



Coherent Behavior, Quasi- Symmetries 
and Quasi-Labels 

Recently the notion of a quasi-symmetry and adiabatic mixing has been introduced in 
nuclear physics [51]. The toy model of a harmonic oscillator in a one-dimensional box can be 
used to introduce and illustrate one possible definition of a quasi-symmetry, an asymptotic 
label (quasi-label), and a coherent behavior associated with a quasi-symmetry. 

First, we define a similarity relation of two states |$) and |\E') with respect to some 
Hermitian operator "H, and denote it as: 

1$) ~ 1^) . 

In this approach, the operational definition of such a similarity relation uses the eigenvectors 
["H; A) and the eigenvalues A of "H: 

H\H;A)=A \ n;A). 

We would say that |$) and |\&) are 7i similar (|$) ~ |^)), if there is a function /, eventually 
monotonic, that maps the distribution p$ (A) = {{n; A\<^)\ to p^, (A) = A|^)f so that: 

K (A)I'^Ip* (/(A))I'. 

In simple words, this means that the shape of the probability distribution p$ is similar 
to the shape of the probability distribution p^. 

If G is a symmetry group for the operator T-L, then the eigenspace for a given A may 
be degenerate, and any function |$) obtained from by a unitary transformation U E G 
(1$) = U l"^)) will be H similar to Thus, G defines an intrinsic symmetry for the wave 
functions that are similar to Therefore, can be viewed as an "intrinsic state" with 
respect to the symmetry of H. 

In the case when T-L is one of the exact limits of a Hamiltonian H, i.e. H = H + X'^V, 
then one can define an adiabatic mixing of the states ["H; A) due to the interaction V. An 
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asymptotic label (quasi-label) A can be assigned to each eigenvector l^'; A) oi H — H + X~^V 
in the limit A — > oo: 

?^ 1^; A, A ^- oo) = A I*; A, A ^ oo) 

thus, 

1^; A, A) 1?/; A) , when A oo. 

Assigning A by using the natural order of the levels must be done carefully by tracing the 
sign of each level crossing. 

Once the asymptotic label A has been assigned for a state \^), then a coherent behavior 
with respect to an observable can be defined as well. There is a quasi-symmetry T-t for the 
observable O : |^) ^ R, if its value O [|\&; A, A)] does not depend much on the parameter A: 

0[|*;A,A)]«0[|?^;A)]. 

Some common functions for O [|^; A, A)] are related to the expectation values of O : 

A, A ^ oo| O I*; A, A ^ oo) (A| O |A) . 

In general, a coherent behavior with respect to other quantities can be defined as well. For 
example, a relative transition rate from a state |^; J) to a state l^'; J + 2) due to a transition 
operator, say E2, can be defined, say by using B {E2, ^, J), where ^ is the "intrinsic state" 
upon which the band is built. All the states within the band should actually be within the 
class of "H equivalent states. In particular, the asymptotic label A can be used as a band 
label. 

Notice that, as in the toy model studied, it may happen that at finite A the wave function 
1^) has been assigned label A while its components along the space |'H;A) are practically 
missing. Following the results from the toy model, we can define some possible types of 
spectral structures that may exhibit such coherent behavior with respect to a Hamiltonian 
H — H + V, and thus to specify a quasi-symmetry. 

Specifically, setting X — 1, H — H + V has a quasi-symmetry if: 

i/|*;A„) = £;(A„)|*;A„), 

'H |A„) = A„ |A„) , A„+i > A„, 

An > {An \V\ An) > A„+i - A„, 
E{An) ^ {An\H\An). 

Here, A„ is the corresponding asymptotic label of the state |^; A„) . The term (A„ \ V\ An) > 
An+i — An means that V mixes strongly different |A„) states. Therefore, perturbation theory 
cannot be applied in the usual small perturbation regime. However, A„ > (A„ \ V\ An) together 
with E (An) ~ (A„ \H\ An) = A„ + (A„ \V\ An) means that the spectral structure of H in this 
region is similar to the spectral structure of "H within a few percent: 

E (An) - An _ (A„|V-|A„) _ {An\V\An) 

E{An) ~ A„ + (AJ\/|A„) ~ An 

This seems to be the situation discussed in the toy model case. Since A„ ~ n^, then 
2A„ > An+i gives {n — 1)^ > 8 which is always satisfied for n > A. 



102 



Appendix E 

Guide to the Oblique-Basis Package 



Although most of the routineq^ in the Oblique-Basis Package 2002 can be compiled by 
using the makefile routines provided, there is a need to follow a few simple steps in order 
to be able to carry out oblique-basis calculations. Here we describe some of the technical 
problems and their solutions that one may face in using the Oblique-Basis package 2002. 
The process of running an oblique-basis calculation consists of four main steps which are 
described below: 

(1) selecting and generating basis states [nuke, PNGGMJ), 

(2) preparing interaction file(s) {IsoInt2pn, M akelnter actions) , 

(3) evaluating matrix elements for the chosen interactions {suSpn), 

(4) solving the generalized eigenvalue problem which includes: obtaining the eigenvectors 
of the Hamiltonian, calculating expectation values and transition probabilities for some 
desired operators {GLanczos). 

In the parentheses are given the names of some relevant routines. Some current limi- 
tations of the Oblique-Basis package 2002 are related to the single-particle basis currently 
employed in the computations. In its present form, the main routine suSpn, which is used to 
evaluate the matrix elements of the operators (step 3 above), is set to operate on spherical 
and cylindrical single-particle states. This clearly restricts the basis generation process (step 
1 above) to the same type of states (spherical and cylindrical). Even though the spherical 
and cylindrical single-particle states are of special interest, the code could, at least in prin- 
ciple, be changed to operate on other desirable single-particle states, such as those one can 
obtain via Hartree-Fock procedure. 

^Thc Oblique Basis package 2002 is available from the author upon request. 
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E.l Generating Basis States 



Since the oblique-basis idea is to combine two or more basis sets, the obhque-code package 
has been designed to use basis states generated and used by other shell-model packages. The 
two main codes in mind are: a version of the Glasgow code, which performs calculations in 
the spherical shell-model basis, and the SU(3) RME code, which performs calculations in 
the SU(3) shell-model basis using cartesian (cylindrical) single-particle states. 

Spherical Shell-Model States. To generate spherical shell-model basis states one needs 
to run the Glasgow code {nuke located in the folder Glasgow — code) with the desired 
configuration hmits and with I BASIS = —1 in the appropriate interaction file (*.mi). For 

more details see the example file BasisOnly.int and files Glasgow.int — instructions or 
Instructions Glasgow test9.doc. There is a script file, RunGlasgow.se, which one may find 
useful when running nuke. This script file uses head files (*.head) and an interaction file 
(*.int) to construct different input files (*.inp) for nuke. 

The wave functions from nuke are stored in the file fort.60 (see Basis — states.info for 
its structure). This file is used by GlsgwBasis2Redstick, located in the folder Oblique — 
Glasgow, to produce file fort.35 which contains the spherical shell-model basis states and 
the single-particle states data. If nuke (Glasgow) is used to calculate the eigenvectors of a 
Hamiltonian, then these eigenvectors are extracted from file fort.60 into file fort.36. 

Important note: File fort.ll is very important pre-generated file. This file (fort.H) 
is used by nuke and GlsgwBasis2Redstick. Do not lose the file fort.lll 

Tip: In oblique runs of the suSpn code, consider using sorted spherical shell-model states. 
Such states are produced by the EpsSorting code located in the folder Oblique/ SphToGyl. 

SU(3) Shell-Model States. Even though the C. Bahri's SU(3) RME code may provide 
the SU(3) highest weight states in the near future, the Oblique -Basis package 2002 has its 
own highest-weight state generator for some of the most important SU(3) irreps. The SU(3) 
related fortran codes are in the folder SUSGenerator. 

The current highest weight state generator SU3JiWS.GEN is located in a folder with 
the same name. Before using SU3.HWS.GEN, one may need to run the auxihary code 
SUSLister which would help in the HWS selection process. The SUSLister code also 
generates the necessary cylindrical single-particle states (file cylin.sps) that are needed in 
the evaluation of matrix elements of operators through the su3pn code. 

Important note: Before doing any runs, compile file fort. A using SU3GENBK located 
in the folder ProjectLibrary . 

Tip: By comparing the output of the SUSLister with the output of the genwsirl, one 
can determine the irreps that are not generated by the generator SUS-HWSjGEN . 

Once the highest weight states are created, they are stored in a file with extension 
*.hws. This file is used as input file to the SU(3) code PNGGMJ located in the sub-folder 
PNGGMJ. The code PNGGMJ generates two files containing basis states *.suS and *.bas. 
It also generates a file PNGGMJ _Brief Jnfo.log which contains information that may be 
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used to set the parameters max_nbas_su3 and mELxmpc for the suSpn code. The *.su3 
file has some of the SU(3) related information and is usually used in testing tools which are 
located in the sub-folder ProjectTools. The *.bas file is mainly for use by the suSpn code. 

E.2 Generating Interaction Files 

Even though there are a few the realistic interactions commonly used (KB3, Wildenthal...), 
their format files may differ significantly. Interactions given in the isospin format as used 
by the Glasgow code can be transformed into the proton-neutron format by using the code 
IsoInt2pn located in the folder Oblique — Glasgow. 

Often used schematic interactions are also available through a package called Makelnter actions 
made by Dr. C. Johnson. M akelnter actions allows one to generated combinations of fre- 
quently studied nuclear interactions. The menu of the currently available interactions is: 

(0) Random noise (TBRE, two-body random ensemble), 

(1) Pairing, 

(2) Multipole-multipole (you choose L), 

(3) S"2 (total spin), 

(4) L"2 (total orbital angular momentum), 

(5) J" 2 (total angular momentum), 

(6) L*S (spin-orbit) 1+2 body. 

E.3 Running the Main Routines 

Evaluating Matrix Elements of Operators. Usually, the running of programs takes 
considerable time. Thus, it is better to use a script file for such runs. Some example script 
files are RunSuSpn + GLDriver.sc and RunSu3pnGLanzos.sc. Basically, the input of the 
suSpn code requires the following entries to be specified: 

> name of the file containing single-particle levels with extension *.sps, 

> scaling of the two-body matrix elements {A/B) " X, 

> interaction file (*.int), 

> name of the file containing cylindrical single-particle levels (*.sps), 

> name of the file containing the su3 basis states (*.bas) 
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> name of the file containing the ssm basis states (*.bas) 

> desired name for the output files (*.ham and *.ovr) 

Eigenvectors, Expectation Values and Transition Matrix Elements. After gener- 
ation of the Hamiltonian and operator matrices, and the overlap matrices (if needed), which 
are stored in files with extensions *.ham and *.ovr, one has to run the generalized Lanczos 
code {GLanczos) to obtain the eigenvalues, eigenvectors, expectation values, and transition 
matrix elements. There is a script file RunGLDriver.se which may be used. 
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